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behavior  extending  over  a  longer  period  for  the  high  strength 
materials.  This  well-defined  time  dependence  clearly  shows 
that  lifetime  predictions  based  on  accelerated  short  hold  time 
tests  could  be  very  non-conservative.  The  influence  of 
humidity  was  also  determined  using  vacuum  as  a  reference 
environment. 

The  second "aspect  of  the  research  investigated  the  nature 
of  the  residual  displacement,  d,  that  leads  to  the  crack  clo¬ 
sure  phenomena.  In  this  aspect  of  the  study,  d  was  determined 
using  a  compliance  gage  on  the  side  of  the  specimen  and  the 
crack  tip  was  monitored  simultaneously  with  an  acoustic  trans¬ 
ducer.  Results  for  varying  AK  in  a  humid  environment  showed 
the  possiblility  of  measuring  d,  using  it  in  a  finite  element 
model  to  predict  the  crack  closure  load.  This  could  then  be 
used  to  predict  AKeff  and  ultimately  fatigue  crack  growth. 

This  trend  was  not  followed  when  the  d  value  was  determined 
in  a  variety  of  environments.  In  fact  an  inverse  correlation 
between  d  and  crack  closure  and  fatigue  crack  growth  rate  v/as 
observed.  This  was  interpreted  in  terms  of  the  unconstrained 
plane  stress  surface  region.  Measurements  are  also  reported 
showing  that  there  is  closure  in  the  interior  plane  strain 
region  that  is  modified  if  the  outer  region  is  machined  away. 
The  results  in  the  d  aspect  of  the  study  was  to  reduce  the 
likelihood  of  the  use  of  surface  measured  d  as  a  useful  set  of 
data  to  use  in  analytic  lifetime  predictions. 

A  The  last  major  emphasis  was  to'  investigate  how  gaseous 
environments  interact  with  fatigue  crack  growth  in  aluminum 
alloys.  Combined  AES/SIMS  studies  were  made  of  fracture  sur¬ 
faces  taken  from  fatigue  crack  growth  specimens  fatigued  in 
vacuum  and  subsequently  exposed  to  and  fatigued  in  a  range  of 
environments.^  The  environments  included  03,  0^8,  H2,  D2,  D2O, 
H2O  and  N2.  r?hq  results  showed  a  deeper  penetration -of  oxygen 
during  fatigue.  \The  deuterium  was  trapped  near  the  surface 
for  the  sample  fatigued  in  D2O  but  not  for  the  sample  fatigued 
in  vacuum  and  subsequently  exposed  to  D2O.  To  help  interpret 
the  results,  deuterium  was  ion  implanted  into  an  aluminum 
specimen  with  varied  surface  conditions  including  the  fatigue 
fracture  surface.  The  diffusion  of  the  deuterium  from  the 
implanted  distribution  was  monitored  by  combined  inert  ion 
sputtering  and  SIMS.  The  deuterium  profiles  obtained  in  the 
fatigue  studies  were  reproduced  by  the  ion  implantation  exper¬ 
iments  leading  to  the  conclusion  that  many  strong  trapping 
sites  were  produced  near  the  fracture  surface  during  the 
fatigue  crack  growth.  A  model  using  an  image  force  technique 
was  developed  to  explain  some  of  the  observed  diffusion 
results  for  the  ion  implanted  deuterium  experiments. 
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THE  INFLUENCE  OF  HOLD  TIMES  ON 
FATIGUE  CRACK  GROWTH  OF  ALUMINUM  ALLOYS 

FINAL  REPORT 
H.L.  Marcus 

Departments  of  Mechanical  Engineering  and 
Materials  Science  and  Engineering 
The  University  of  Texas 
Austin,  Texas  78712 

I .  Introduction 

This  report  will  describe  the  results  of  a  series  of 
experiments  whose  primary  aim  was  to  clarify  the  fatigue 
crack  growth  behavior  of  aluminum  alloys.  The  results  are 
presented  in  a  condensed  form  in  the  main  body  of  the  report 
with  the  details  of  the  results  in  a  series  of  appendices, 
which  in  general  are  either  preprints  or  reprints  of  papers 
submitted  for  publication.  The  major  literature  references 
are  also  included  in  the  write-ups  in  the  appendices. 

There  were  several  aspects  of  fatigue  crack  growth  in 
aluminum  alloys  that  were  investigated.  The  first  was  the 
study  of  the  influence  of  hold  time  at  an  underload  after 
an  overload  on  the  crack  retardation  behavior  of  the  aluminum 
alloys.  As  part  of  this  effort,  the  crack  closure  behavior 
of  the  material  was  observed  using  ultrasonic  and  compliance 
techniques.  To  clarify  the  ultrasonic  technique,  a  study  on 
the  acoustic  energy  transmissions  through  an  interface  was 

made.  A  second  part  of  the  study  was  an  attempt  to  identify 
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the  cause  of  crack  closure  in  terms  of  the  nature  of  the 
material  in  actual  contact.  An  added  incentive  to  getting  a 
direct  measure  of  the  origin  of  crack  closure  by  a  measure¬ 
ment  of  the  amount  of  residual  displacement  was  its  potential 
use  as  a  characteristic  number  for  closure  load  prediction. 
This  could  be  used  for  crack  growth  prediction  using  finite 
element  methods  to  define  AK^f . 

In  order  to  evaluate  this  behavior,  a  series  of  experi¬ 
ments  were  performed  in  hard  vacuum  and  a  range  of  gaseous 
environments.  The  aluminum  alloys  2219-T851  and  7075-T651 
were  used  in  the  studies  in  the  above  conditions  and  in  the 
lower  strength  level  overaged  condition. 

To  further  characterize  the  fracture  surfaces  a  series 

of  secondary  ion  mass  spectroscopy  (SIMS)  and  scanning 

auger  microscopy  (SAM)  experiments  were  performed  on  samples 

18 

fatigued  in  gaseous  environments  including  0  ,  D2°'  D2  aru* 

vacuum.  Additional  ion  implantation  SIMS  studies  were  done 
to  assist  in  evaluating  the  observed  fatigue  results. 

II .  The  Influence  of  Underload  Time  on  Overload  Crack 

Retardation 

It  is  well  known  that  an  overload  during  a  constant 
load  amplitude  fatigue  crack  growth  study  will  greatly 
retard  the  subsequent  crack  growth.  In  addition,  an  under¬ 
load  following  the  crack  growth  reduces  the  magnitude  of 
the  retardation  phenomenon.  In  this  study  the  influence  of 
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hold  time  at  an  underload  after  an  overload  was  investigated 
in  7075-T651  and  2219-T851  and  in  their  overaged  condition 
for  a  range  of  loading  and  environmental  conditions.  The 
details  of  these  results  are  presented  in  Appendices  A  and  B. 
The  main  results  of  the  underload  study  were  the  following: 

1)  A  significant  reduction  in  retardation  due  to  an  overload 
is  observed  for  a  zero  load  underload  when  the  underload  is 
held  for  24  to  48  hours.  If  this  hold  time  effect  is  not  taken 
into  account,  the  predicted  lifetime  based  on  the  accelerated 
short  underload  time  tests  can  be  more  than  twice  the  actual 
life. 

2)  The  time  required  for  the  underload  to  complete  its  relax¬ 
ation  is  a  function  of  the  strength  level  of  the  material. 

For  both  the  2219-T851  and  the  7075-T651  alloys  the  relaxation 
effect  due  to  the  underload  is  approximately  complete  after  24 
hours  but  is  significantly  less  when  the  material  is  overaged 
to  significantly  reduce  the  yield  strength. 

3)  There  is  a  difference  in  the  underload  hold  time  effect 
depending  on  whether  the  experiment  was  run  in  vacuum  or  in 
humid  air  (relative  humidity  between  5  and  50%)  with  a 
greater  relaxation  effect  seen  in  humid  air. 

4)  The  retardation  seems  to  be  related  to  crack  closure 
resulting  from  the  crack  tip  being  enclosed  by  the  residual 
deformation  and  displacement  from  the  overload.  Hold  time  at 
underloads  apparently  reduces  the  magnitude  of  the  residual 
displacements  through  some  room  temperature  relaxation  effect. 


-4- 


III .  Ultrasonic  Transmission  Through  A  Rough  Interface 

In  the  studies  reported  here  the  measurement  of  crack 
closure  using  ultrasonic  transducers  is  heavily  emphasized. 

One  of  the  major  questions  about  this  technique  to  identify 
the  closure  load  is  how  much  force  across  the  fracture  surface 
is  required  to  transmit  the  acoustic  wave.  In  principle,  the 
transmission  across  a  very  thin  air  gap  can  be  significant. 

To  determine  the  stresses  required,  a  series  of  surfaces  of 
the  aluminum  alloys  were  prepared  with  a  wide  range  in  sur¬ 
face  roughness.  The  ultrasonic  signal  transmitted  through 
this  interface  was  then  measured  as  a  function  of  applied 
load.  The  results  clearly  demonstrated  a  normal  stress  of 
the  order  of  five  to  twenty  percent  of  the  yield  strength  was 
required  for  significant  transmission  of  the  ultrasonic  wave. 
Virtually  no  signal  was  transmitted  when  the  blocks  lay  at 
rest  on  each  other.  Based  on  these  results,  the  details  of 
which  are  given  in  appendix  C,  the  ultrasonic  closure  measure¬ 
ments  were  considered  to  be  realistic  since  a  significant  com¬ 
pressive  stress  was  required  for  transmission  of  the  signal. 

No  attempt  was  made  to  evaluate  diffraction  effects  on  the 
transmitted  ultrasonic  wave  associated  with  the  change  in 
crack  tip  shape  during  the  loading  cycle. 

As  an  outgrowth  of  these  measurements  a  model  was  developed 
to  identify  the  true  crack  length  using  ultrasonic  measurements 
when  the  crack  is  partially  closed.  To  accomplish  this 


nondestructive  evaluation  of  the  crack  length,  the  stress  and 
roughness  dependence  of  the  transmitted  acoustic  energy  is 
integrated  to  give  the  measured  ultrasonic  intensity.  This 
approach  could  be  used  in  either  the  transmitted  or  back- 
scattered  mode.  The  details  of  the  model  are  given  in 
appendix  D  where  several  specific  crack  and  transducer  geo¬ 
metries  are  discussed. 

IV.  Evaluation  of  the  Residual  Displacement,  d.  Associated 

With  Crack  Closure 

In  the  Elber  model  of  crack  closure  the  two  fracture 
surfaces  can  come  in  contact  during  the  low  end  of  the  ten¬ 
sion  cycle  only  if  displacement  above  the  theoretical  plane 
of  the  crack  exists.  This  residual  displacement,  d,  has 
several  potential  origins  including  surface  roughness,  volume 
increase  due  to  deformation  including  the  constant  volume 
effect  and  formation  of  an  oxide  on  the  surface.  The  con¬ 
stant  volume  effect  associated  with  the  plastic  deformation 
will  be  most  significant  in  the  vicinity  where  the  crack  tip 
intercepts  the  surface  of  the  sample.  In  this  region  the 
plane  strain  condition  is  relaxed  and  the  behavior  is  more 
closely  simulated  by  a  plane  stress  description. 

In  considering  how  to  analyze  the  residual  displacement, 
several  approaches  were  taken.  In  the  first  approach,  it  was 
assumed  that  if  you  could  measure  the  d  value,  it  would  be 
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representative  of  the  crack  closure  behavior.  If  this  was 
true,  the  crack  closure  load  could  be  calculated  using  a  finite 
element  model  where  the  measured  residual  displacement  was 
introduced  onto  the  idealized  fracture  surface.  A  description 
of  the  measurement  of  d  using  surface  clip-on  gauges  and  the 
finite  element  model  is  given  in  appendix  E.  In  this  case  the 
value  of  d  was  determined  in  a  humid  environment  and  the  d 
value  determined  was  found  to  increase  with  the  closure  load. 
Good  agreement  was  found  over  a  limited  AK  range  between  the 
ultrasonically  measured  crack  closure  load  and  the  crack  clo¬ 
sure  load  calculated  from  the  finite  element  model. 

If  this  correlation  between  d  and  the  crack  closure  load 
could  be  confirmed  over  a  wide  range  of  loading  and  environ¬ 
mental  conditions,  the  possibility  would  then  exist  to  measure 
d  for  any  expected  load  history,  and  using  the  finite  element 
analysis  predict  the  closure  load  and  directly  be  able  to  cal¬ 
culate  life.  In  order  to  more  closely  examine  this  correlation 
a  series  of  measurements  were  performed  in  vacuum  and  in  humid 
air.  The  details  of  these  experiments  are  given  in  appendix  F. 
There  were  several  conclusions  of  interest  coming  from  these 
results.  The  first  was  the  fact  that  the  measured  value  of  d 
determined  by  running  the  crack  in  the  compact  tension  specimen 
through  a  clip-on  gage  mounted  on  the  surface  was  inversely 
proportional  to  the  crack  closure  load  measured  ultrasonically 
and  the  crack  growth  rate  as  a  function  of  relative  humidity. 


This  was  totally  opposite  to  what  was  expected.  To  evaluate 
this,  a  series  of  experiments  were  run  measuring  d  and  crack 
closure  where  the  crack  was  grown  in  vacuum  and  then  various 
gaseous  environments  were  introduced.  It  was  found  that  the 
d  value  increased  dramatically  when  water  vapor  was  introduced 
but  the  ultrasonically  measured  crack  closure  load  was  lowered. 
The  introduction  of  N2  had  no  effect  and  02  a  very  small  effect. 
The  results  led  to  the  conclusion  that  the  humidity  influenced 
the  near  surface  plane  stress  region  to  a  much  greater  extent 
than  the  interior.  A  possible  explanation  is  the  oxide 
strengthening  of  the  asperities  led  to  the  increase  in  d. 

This  sequence  of  results  clearly  demonstrated  that  the 
concept  of  using  surface  measured  d  values  to  predict  overall 
crack  closure  was  not  a  viable  approach.  In  fact  it  pointed 
out  the  extreme  caution  that  must  be  used  in  interpreting  any 
surface  strain  measurement. 

Another  factor  in  the  measurement  of  crack  closure  is  the 
question  of  whether  it  is  only  a  plane  stress  phenomena.  To 
clarify  this  point,  crack  closure  was  measured  on  a  sample 
using  ultrasonic  methods.  The  outsides  of  the  sample  were 
then  machined  away  and  the  closure  load  remeasured.  The 
details  of  these  measurements  are  given  in  appendix  G.  Crack 
closure  was  present  in  the  ultrasonic  measurements  in  the  cen¬ 
tral  plane  strain  region  of  the  specimen  before  machining. 

The  crack  closure  load  was  decreased  when  the  outer  plane 


stress  region  with  the  high  d  was  removed.  This  was  due  to 
the  lack  of  the  high  d  outer  region  propping  open  the  center 
of  the  crack  during  the  unloading  cycle.  The  results  clearly 
demonstrate  that  crack  closure  occurs  in  the  plane  strain 
region. 

V.  The  Interaction  of  Gases  with  the  Crack  in  Aluminum  During 
Fatigue  Crack  Growth 

The  final  area  of  investigation  probing  the  fatigue 

crack  growth  in  aluminum  alloys  was  the  attempt  to  evaluate 

the  interaction  between  the  gaseous  environments  and  the 

fatigue  crack.  In  this  study  an  attempt  was  made  to  identify 

where  the  elements  went  relative  to  the  surface.  The  fatigue 

experiments  consisted  of  running  a  fatigue  crack  in  vacuum 

for  a  short  way  and  then  introducing  the  gaseous  environment 

and  continuing  the  fatigue  crack  growth.  The  environments 
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studied  included  vacuum,  C^,  ^ ,  I^O,  D20,  °2 '  °2  '  an^  °2  * 

The  samples  were  removed  from  the  mechanical  test  chamber 
after  fatigue  crack  growth  was  completed  and  transferred  to 
the  SAM/SIMS  system.  The  details  for  these  studies  are  given 
in  appendices  H,  I  and  J.  In  all  cases  inert  ion  sputtering 
combined  with  AES  and  SIMS  was  performed  in  the  fracture  sur¬ 
face  formed  in  vacuum  and  subsequently  exposed  to  the  gas  and 
on  the  surface  formed  in  the  gas.  In  the  case  of  oxygen,  it 
was  consistently  found  that  the  oxygen  penetrated  deeper  into 
the  fatigue  crack  fracture  surface  than  the  normal  oxide 


thickness  associated  with  diffusion  (appendix  H) .  This  was 
true  for  the  oxide  formed  in  both  gaseous  oxygen  or  in  the 
presence  of  water  vapor.  This  extra  penetration  from  dry 
oxygen  did  not  seem  to  enhance  the  fatigue  crack  growth  rate 
in  the  aluminum  alloys,  although  it  did  in  the  Ti  and  Ni/Cu 
alloys  that  were  also  tested.  This  lack  of  enhanced  fatigue 
crack  growth  in  the  aluminum  alloys  exposed  to  dry  oxygen 
even  though  the  oxygen  was  deeper  into  the  surface  could  be 
an  indication  that  in  th°  case  of  water  vapor  the  enhanced 
crack  growth  rate  is  probably  not  associated  with  the  enhanced 
oxygen  presence.  This  would  leave  the  hydrogen  originating 
from  the  water  as  the  element  most  likely  enhancing  the 
fatigue  crack  growth  of  the  aluminum  alloys. 

In  order  to  evaluate  the  region  of  hydrogen  interaction, 
SAM/SIMS  studies  were  performed  on  the  fracture  surface  of 
the  alloys  fatigued  first  in  vacuum  and  then  in  D2O.  The 
deuterium  was  used  as  a  marker  to  control  hydrogen  contami¬ 
nation  arising  during  the  transfer  between  the  fatigue  envi¬ 
ronmental  chamber  and  the  SAM/SIMS  vacuum  chamber.  The  SAM 
was  used  to  monitor  the  oxygen  peak  and  the  SIMS  the  deuterium 
peak  during  the  inert  ion  sputtering  (appendix  I) .  The 
results  for  the  D2O  measurements  show  a  much  greater  penetra¬ 
tion  of  deuterium  into  the  fracture  surface  created  during 
exposure  to  D2O  than  the  segment  exposed  to  D2O  after  the 
fatigue  crack  was  formed. 
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The  results  were  obtained  about  24  hours  after  the 
fatigue  crack  growth  took  place.  If  there  were  no  strong 
trapping  sites  available  for  the  deuterium,  it  should  have 
diffused  out  of  the  aluminum  alloys  during  the  time  of 
fatiguing  and  transfer.  In  order  to  investigate  the  pres¬ 
ence  of  trapping  sites,  a  series  of  deuterium  ion  implanta¬ 
tion  experiments  were  performed  on  the  7075  and  2219  alloys 
that  were  initially  in  the  unrolled,  rolled,  fatigue  frac¬ 
ture  surface  and  unrolled  but  damaged  by  high  energy  neon 
ion  implantation  prior  to  the  time  of  the  deuterium  ion 
implantation.  For  the  details  of  the  results  and  analysis, 
see  appendix  J.  The  deuterium  profiles  obtained  could  be 
broken  into  two  categories.  The  first  was  for  the  unrolled 
specimen,  with  and  without  neon  ion  induced  damage.  In 
these  cases  the  results  were  controlled  by  the  surface  oxide 
acting  as  a  diffusion  barrier.  To  interpret  this  behavior, 
an  image  force  model  for  diffusion  across  the  oxide  surface 
was  developed.  The  results  of  the  model  closely  correlated 
with  the  experimental  results.  The  second  type  of  result  was 
obtained  for  the  fatigued  fracture  surfaces  and  the  cold 
worked  surfaces.  In  this  case,  the  deuterium  profiles  were 
very  similar  in  depth  beyond  the  oxide  to  that  found  in  the 
samples  fatigued  in  D2O.  The  profiles  in  the  fatigued  frac¬ 
ture  surfaces  associated  with  the  implanted  deuterium  were 
formed  by  diffusion  toward  the  surface  with  the  trapping  sites 
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filling  during  the  exit  of  the  deuterium.  The  samples 
fatigued  in  D2O  resulted  in  deuterium  profiles  where  the 
deuterium  originates  at  the  crack  tip.  The  fact  that  both 
resulted  in  similar  profiles  can  be  interpreted  to  mean  that 
the  observed  distribution  is  representative  of  the  strong 
trapping  sites  distribution.  This  distribution  during 
fatigue  can  influence  the  fatigue  crack  growth  rate,  but 
the  mechanism  of  its  influence  is  not  modeled  at  present. 

VI .  Summary 

The  results  of  the  research  on  fatigue  crack  growth  in 
aluminum  alloys  are  described  in  three  general  areas.  The 
first  relates  the  influence  of  hold  time  at  an  underload  that 
follows  an  overload  on  the  fatigue  crack  growth  retardation 
behavior  of  the  7075  and  2219  aluminum  alloys.  A  very  dis¬ 
tinct  effect  on  retardation  with  hold  time  is  observed.  The 
effect  is  microstructure  dependent  with  the  time  dependence 
behavior  extending  over  a  longer  period  for  the  high  strength 
materials.  This  well-defined  time  dependence  clearly  shows 
that  lifetime  predictions  based  on  accelerated  short  hold  time 
tests  could  be  very  non-conservative.  The  influence  of 
humidity  was  also  determined  using  vacuum  as  a  reference 
environment. 

The  second  aspect  of  the  research  investigated  the  nature 
of  the  residual  displacement,  d,  that  leads  to  the  crack  clo¬ 
sure  phenomena.  In  this  aspect  of  the  study,  d  was  determined 
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using  a  compliance  gage  on  the  side  of  the  specimen  and  the 
crack  tip  was  monitored  simultaneously  with  an  acoustic  trans¬ 
ducer.  Results  for  varying  AK  in  a  humid  environment  showed 
the  possibility  of  measuring  d,  using  it  in  a  finite  element 
model  to  predict  the  crack  closure  load.  This  could  then  be 
used  to  predict  and  ultimately  fatigue  crack  growth. 

This  trend  was  not  followed  when  the  d  value  was  determined 
in  a  variety  of  environments.  In  fact  an  inverse  correlation 
between  d  and  crack  closure  and  fatigue  crack  growth  rate 
observed.  This  was  interpreted  in  terms  of  the  unconstrained 
plane  stress  surface  region.  Measurements  are  also  reported 
shoving  that  there  is  closure  in  the  interior  plane  strain 
region  that  is  modified  if  the  outer  region  is  machined 
away.  The  results  in  the  d  aspect  of  the  study  was  to 

reduce  the  likelihood  of  the  use  of  surface  measured  d  as  a 
useful  set  of  data  to  use  in  analytic  lifetime  predictions. 

The  last  major  emphasis  was  to  investigate  how  gaseous 
environments  interact  with  fatigue  crack  growth  in  aluminum 
alloys.  Combined  AES/SIMS  studies  were  made  of  fracture  sur¬ 
faces  taken  from  fatigue  crack  growth  specimens  fatigued  in 

vacuum  and  subsequently  exposed  to  and  fatigued  in  a  range  of 

18 

environments.  The  environments  included  C^,  C>2  ,  D2,  D20' 

H20  and  N2 •  The  results  showed  a  deeper  penetration  of  oxygen 
during  fatigue.  The  deuterium  was  trapped  near  the  surface 
for  the  sample  fatigued  in  D2O  but  not  for  the  sample  fatigued 


in  vacuum  and  subsequently  exposed  to  D2O.  To  help  inter¬ 
pret  the  results,  deuterium  was  ion  implanted  into  an  alumi¬ 
num  specimen  with  varied  surface  conditions  including  the 
fatigue  fracture  surface.  The  diffusion  of  the  deuterium  from 
the  implanted  distribution  was  monitored  by  combined  inert 
ion  sputtering  and  SIMS.  The  deuterium  profiles  obtained  in 
the  fatigue  studies  were  reproduced  by  the  ion  implantation 
experiments  leading  to  the  conclusion  that  many  strong  trapping 
sites  were  produced  near  the  fracture  surface  during  the 
fatigue  crack  growth.  A  model  using  an  image  force  technique 
was  developed  to  explain  some  of  the  observed  diffusion 
results  for  the  ion  implanted  deuterium  experiments. 
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It  is  now  well  known  that  an  overload  applied  during  a  metal  fa¬ 
tigue  crack  growth  experiment  causes  crack  growth  retardation  (1,2,31. 

In  addition,  underloads  subsequent  to  the  overload  tend  to  reduce  the 
retardation  phenomena  [1].  It  has  been  reported  that  time  at  small 
underload  does  not  seem  to  play  a  significant  role  in  its  influence  on 
the  retardation  behavior  [3].  The  data  presented  here,  however,  where 
t he  underload  is  greater,  do  show  such  a  time  effect.  This  question  is 
of  particular  significance  in  terms  of  relating  accelerated  laboratory 
tests  to  the  actual  loading  seen  during  the  service  lifetime  of  a  struc¬ 
ture  designed  to  the  fatigue  crack  growth  behavior  of  the  material. 


Fatigue  crack  growth  experiments  were  carried  out  on  an  A1  2219- 
T851  (oy  =  50  ksi)  compact  tension  specimen  ( B  =  1  in,  W  =  4.25  in, 

H  =  2.09  in).  Baseline  tension-tension  loading  consisted  of  a  10  Hz 
sine  wave  whose  mean  value  and  amplitude  are  varied  so  that  at  each 
point  where  an  overload  is  applied  (4  points  per  sample,  spaced  3/8  in 
apart),  Kmjn  =  4  ksi/in  and  Kmax  *  12  ksi/in.  The  overload  is  such 

that  KOVOri0ad  =  30  ksi,/in  (K|  2  40  ksi/in).  Testing  is  done  on  an 

MTS  closed  loop  electro-hydraul ic  testing  machine,  in  room  air  with 
relative  humidity  between  30%  and  55%. 


Crack  depth  is  plotted  versus  number  of  cycles.  The  extension  of 
the  steady-state  growth  lines  before  and  after  the  overload  should  be 
roughly  parallel,  with  a  separation  of  a  certain  number  of  cycles  which 
we  will  call  the  delay. 


The  relative  crack  depth  is  measured  by  the  attenuation  of  a  trans¬ 
mitted  longitudinal  ultrasonic  wave,  transmitted  and  received  by  two 
2.25  MHz  piezoelectric  transducers.  The  received  signal  is  converted 
to  a  DC  signal  and  is  plotted  against  load.  This  plot  is  not  only  a 
means  of  monitoring  crack  growth,  but  also  shows  the  crack  closure  ef¬ 
fect  reported  in  (2,3,4),  that  is,  that  the  crack  effectively  closes  at 
loads  below  some  l’c  >  Pmjn. 

An  underload  of  varying  magnitude  and  duration  was  applied  after 
each  overload.  Each  test  was  repeated  three  times,  never  twice  on  the 
same  sample  or  at  the  same  crack  length.  Results  are  summarized  in 
Figs.  I  and  2. 


The  extent  of  the  retardation  delay  decreases  with  the  time  at  zero 
load  (Fig.  1)  and  also  with  the  magnitude  of  the  underload  (Fig.  2). 

The  significant  result  shown  here  is  the  influence  of  holding  time  at 
zero  load.  A  decrease  of  about  50%  in  retardation  delay  is  caused  by 
an  underload  of  about  2  sec.  Extending  the  hold  time  to  24  hr  reduces 
this  by  about  another  30%.  The  impact  on  design  philosophy  would  be 
even  greater  if  the  2  sec  hold  time  represents  the  accelerated  tests 
and  the  24  hr  represents  the  real  loading  history. 
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tKing  the  crack  closure  model  [2,4],  what  happens  is  that  when  base¬ 
line  loading  is  resumed,  the  crack  is  not  fully  open  and  the  crack 
growth  rate  is  reduced  and  will  not  return  to  the  steady-state  value  un¬ 
til  the  two  rough  crack  surfaces  created  during  the  overload  become 
smooth  through  cyclic  contact.  During  the  underload,  these  two  surfaces 
are  smoothed  as  they  are  pressed  together.  Time  at  underload  then  in¬ 
creases  this  smoothing  due  to  a  strain  relaxation  process. 

Additional  experimental  effort  is  now  being  performed  to  determine 
the  influence  of  strength  level,  microstructure,  humidity,  and  other 
loading  histories  on  this  relaxation  behavior,  and  an  attempt  will  be 
made  to  model  the  strain  relaxation  phenomena. 
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Introduction 

The  phenomenon  of  crack  retardation  following  an  over- 

f  1-4  1 

load  has  been  studied  extensively  .  It  is  also  reported 

in  the  literature  that  an  underload  following  an  overload 
tends  to  reduce  the  retardation  phenomenon^4  Effects  of 

hold  times  on  this  reduction  in  delay  have  also  been  investi- 
gated  to  a  certain  extent1  .  Such  time  effects  could  be 
significant  in  correlating  laboratory  data  from  accelerated 
life  tests  to  real  life  fatigue  performance  of  engineering 
structures. 

In  this  investigation,  hold  time  effects  on  fatigue 
crack  retardation  in  two  high  strength  aluminum  alloys,  2219 
and  7075,  were  studied.  Influences  of  different  metallurgical 
conditions  and  different  environments  were  examined. 

Many  of  the  observations  in  this  research  will  be  dis- 

r  gi 

cussed  in  terms  of  crack  closure  model  proposed  by  Elber  J 
and  in  terms  of  strain  relaxation  effects  proposed  elsewhere 
[9,10]  ^  ^as  £een  suggested  that  underloads  following  over¬ 


loads  reduce  retardation  due  to  strain  relaxation  at  the  crack 
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tip  .  This  phenomenon  can  take  place  due  to  dislocation 

cross  slip^^.  In  that  case  longer  hold  times  should  cause 

shorter  delay  periods  which  was  observed  by  Sewell .  But 
r  7 1 

Sharpe  et  al  found  only  a  slight  correlation  between  hold 
time  at  zero  load  and  the  number  of  delay  cycles  [ND]  while 
Chanani^^  reported  no  trend  between  and  hold  time  at 
mean  load  following  an  overload. 

An  aggressive  environment  can  affect  the  extent  of  plastic 
deformation.  Retardation  is  expected  to  increase  in  vacuum  or 
dry  atmospheres  as  compared  to  moist  air.  This  trend  was 
reported  by  Sewell  ^4].  Buck  et  al  observed  lower 

values  of  AKgff  in  dry  nitrogen  than  in  moist  air. 

Sewell  and  Marcus ^  report  experiments  where  each  over¬ 
load  was  followed  by  one  hour  relaxation  at  an  underload  of 
variable  magnitude.  The  results  indicated  that  there  was  a 
reduction  in  delay  with  increasing  magnitude  of  underload. 

In  low  humidity  (<5%RH)  experiments,  they  observed  that  all 
the  phases  of  crack  growth  experiments  slowed  down.  These 
include  steady  state  crack  growth  and  crack  growth  following 
any  overload-under load  combination.  A  decrease  in  delay  with 
increase  in  holding  time  at  zero  load  was  again  evident. 

Because  of  these  conflicting  results  it  became  necessary 
to  study  the  effect  of  hold  times  on  the  retardation  phenome¬ 
non.  This  investigation  concerns  with  hold  time  experiments 
in  vacuum  and  humid  air  for  aluminum  alloys  with  different 
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metallurgical  treatments.  Crack  growth  and  retardation  data 
was  obtained  for  various  hold  times.  The  results  are  then 
interpreted  in  terms  of  crack  closure  model. 

Experimental  Procedure 

Experiments  were  performed  on  the  aluminum  alloys  7075- 
T651  and  2219-T851,  as  well  as  overaged  7075  and  2219.  The 
mechanical  properties  are  given  in  Table  1.  Experiments 
were  carried  out  in  humid  air  and  in  a  hard  vacuum  system 
(0.1  pPa) . 

The  experiments  were  performed  on  compact  tension  speci¬ 
mens  with  B  =  2.54  cm.,  W  =  10.8  cm  and  H  =  5.3  cm.  The 
specimens  had  a  T-L  orientation.  An  MTS  closed  loop  electro- 
hydraulic  fatigue  testing  machine  was  used  for  loading.  All 
the  cycling  was  sinusoidal  tension-tension  at  a  frequency  of 
10  Hz  for  a  constant  range  in  stress  intensity,  &K,  and  load 
ratio,  R.  Crack  closure  and  crack  growth  was  measured  using 

2.25  MHz  longitudinal  ultrasonic  waves  in  transmission.  The 

(9  12  131 

details  of  the  technique  are  described  elsewhere  '  '  . 

A  clip  on  gage  at  the  crack  mouth  was  also  used  to  monitor 
closure  and  growth. 

The  actual  relaxation  experiments  were  performed  by  first 
growing  the  crack  for  at  least  48  mm  before  each  overload. 

This  allowed  the  crack  to  reach  its  steady  state  crack  growth 
rate  for  the  applied  AK  and  allowed  the  plastic  zone  to  sta¬ 
bilize.  During  the  preliminary  growth,  the  ultrasonic  intensity 
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was  periodically  recorded  along  with  the  number  of  cycles 
N.  At  a  predetermined  crack  length,  an  overload  KQL,  was 
applied  at  approximately  0.1  Hz  followed  immediately  by  an 
underload  to  zero.  The  load  was  maintained  at  zero  for  a 
specified  hold  time  and  then  raised  to  the  mean  load  where¬ 
upon  baseline  cycling  was  resumed.  Immediately  following 
the  overload,  crack  retardation  was  observed  and  the  ultra¬ 
sonic  signal  would  remain  constant.  After  several  thousand 
cycles  crack  growth  would  start  again  and  the  retardation 
measurements  would  be  made  until  steady  state  da/dN  was 
resumed.  The  number  of  delay  cycles,  ,  was  then  deter¬ 
mined  using  a  tangent  method . 

Four  overloads  spaced  9.5  mm  apart  were  applied  to  each 

specimen.  This  distance  was  chosen  on  the  basis  of  observa- 
.  r  14 1 

tions  by  Wei  J  which  indicates  no  interaction  effect  for  a 
distance  greater  than  13  times  the  size  of  the  plastic  zone. 

To  minimize  any  possible  interaction  effects,  the  tests  were 
varied  so  that  no  hold  time  was  repeated  on  the  same  specimen 
or  at  the  same  crack  length.  The  overloads  were  applied  at 
crack  lengths  of  2.7  cm,  3.65  cm,  4.6  cm  and  5.56  cm  so  that 
all  tests  were  performed  in  the  range  0 . 25<a/w<0 . 51. 

Load  Parameters 

1.  Humid  air:  The  7075  alloy  experiments  were  conducted  with 

loading  parameters  of  K  .  =  3.3  MPa/m,  K _ =  9.9  MPa/m  and 

mm  max 
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Kql  =  22  MPa/m.  The  2219  alloy  experiments  were  performed 

with  K  ■  =  4.9  MPa /in,  K  =13.2  MPa/in  and  K_T  =  29.7  MPa/in. 

min  max  OL 

Hold  times  were  with  K  =  0. 

2.  Vacuum:  The  retardation  experiments  were  repeated  in 

0.1  yPa  vacuum  for  7075-T651  alloy.  The  load  parameters 

were  K  .  =  1.1  MPa/in,  K  =  11  MPa/in  and  KOT  =  22  MPa/in. 

min  max  OL 

This  larger  value  of  AK  was  used  since  crack  growth  was 
significantly  slower  in  the  vacuum  environment. 

Results  and  Discussion 

Humid  air:  The  data  from  7075-T651  alloy  experiments  appears 
in  Fig.  1.  Without  an  underload,  the  delay  times  were  so  long 
(over  300  K  cycles)  that  the  tests  were  aborted  before  the 
crack  started  propagating  again.  These  long  delay  times  were 
probably  due  to  large  plastic  deformation  at  the  crack  tip 
caused  by  the  KQL  which  was  very  close  to  Klc  of  the  material. 
Application  of  an  underload  to  zero  immediately  following  an 
overload  caused  a  great  reduction  in  the  number  of  delay  cycles 
and  holding  the  7075  specimen  at  the  underload  reduced  ND  even 
more.  A  steady  reduction  in  retardation  was  observed  as  the 
hold  time  was  increased,  and  it  seemed  to  level  off  after  24 
hours. 

In  the  overaged  7075,  most  of  the  reduction  in  the  retarda 
tion  due  to  overload  occurred  immediately.  A  two-second  hold 
time  at  zero  load  caused  a  significant  decrease  in  delay.  Only 
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a  small  further  decrease  was  observed  as  the  hold  time 
increased.  The  indication  was  that  relaxation  effects  occurred 
more  quickly  in  the  lower  strength  7075  overaged  material. 

The  data  for  2219  alloys  is  presented  in  Fig.  2.  The 
results  show  that  for  both  2219-T651  and  overaged  2219  alloys, 
a  substantial  decrease  in  retardation  due  to  the  application  of  an 
underload  after  the  spike  is  observed.  This  trend  was  similar  to 
that  in  7075  alloys.  The  ND  values  for  no  underload  for  the 
overaged  2219  were  less  than  those  for  the  T851.  Reduction  in 
Nd  for  the  hold  tine  experiments  showed  a  similar  trend. 

Vacuum:  The  overload  followed  by  zero  load  hold  time  data  for 

7075-T651  alloys  inO.lyPa  vacuum  is  presented  in  Fig.  3.  Again, 
as  in  the  humid  air  tests,  there  was  a  significant  reduction 
in  delay  cycles  with  increasing  time  at  underload. 

More  relaxation  seemed  to  occur  in  humid  air  hold  time 
tests  than  in  vacuum  tests.  This  may  have  been  partially  due 
to  differences  in  loading  conditions  of  vacuum  and  humid  air, 
but  nevertheless  the  results  do  indicate  that  relaxation  effects 
are  significant  in  vacuum  but  are  not  as  great  as  in  humid  air. 

Strain  Relaxation 

The  above  results  can  be  explained  in  terms  of  strain 
relaxation  in  the  vicinity  of  the  crack  tip.  The  dislocation 
cross  slip  mechanism  which  results  in  relaxation,  is  less 
active  in  vacuum  than  in  more  aggressive  environments  like 
humid  air .  Overaging  causes  a  significant  decrease  in  the 
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delay  times  compared  to  the  high  strength  probably  due  to 
the  residual  strain  at  the  crack  tip  relaxing  out  more 
easily  for  the  lower  strength  materials.  Sanders  et  al^*^ 
suggest  that  because  of  coarsening  of  the  dispersoids,  over¬ 
aging  may  reduce  the  number  of  secondary  cracks  at  the  main 
crack  front  resulting  from  fracture  at  particle  matrix  inter¬ 
face.  This  decrease  in  crack  branching  causes  stress  inten¬ 
sity  at  the  main  crack  tip  to  increase  resulting  in  faster 
crack  propagation  and  decrease  in  delay. 

Crack  Closure 

If  strain  relaxation  is  occurring  in  the  vicinity  of 
crack  tip  during  overload-under load  cycle,  then  it  should 
cause  a  noticeable  change  in  crack  closure  load.  To  verify 
this,  experiments  were  carried  out  on  overaged  2219  alloys  with 
load  parameters  mentioned  earlier.  The  ultrasonic  observations 
are  shown  in  Fig.  4.  These  closure  curves  indicate  that  at  the 
beginning  of  the  delay  period  immediately  following  the  overload- 
underload  sequence,  the  closure  curve  is  almost  completely 
vertical  (N  =  171000).  Here  the  crack  tip  is  completely  closed, 
but  the  material  which  has  been  plastically  strained  by  the 
overload  is  preventing  the  crack  surfaces  from  closing  any  fur¬ 
ther.  As  cycling  continues,  the  shape  of  the  curves  begin  to 
change.  The  curves  first  begin  to  shift  to  the  left  indicating 
that  the  tip  is  just  beginning  to  open  at  the  higher  position 
of  the  load  cycle.  Two  inflexion  points  become  apparent 
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indicating  the  crack  closes  a  short  distance.  Soon  the 
bottom  of  the  curve  smooths  out.  A  well  defined  closure  load 
becomes  evident  and  it  decreases  with  cycling  to  pre-overload 
value. 

The  closure  data  presented  here  seems  to  indicate  that 
overload  strained  material  which  has  been  wedging  the  crack 
surfaces  apart  has  been  gradually  hammered  down.  If  this 
noticeable  change  in  closure  curves  does  imply  hammering,  it 
should  show  up  on  the  fracture  surface  also.  To  look  for  such 
hammering,  cross  sections  of the  fracture  surface  of  the  2219- 
T851  specimen  were  examined  with  an  optical  microscope.  The 
overload  region  appears  to  be  smoothed  down  or  smashed  down 
as  compared  to  fast  fracture  region.  This  gives  some  visual 
evidence  of  hammering. 

Thus  delay  after  overload-under load  cycle  is  associated 
with  noticeable  change  in  closure  behavior  of  the  crack  tip. 
Although  it  is  not  clear  whether  the  closure  observations  were 
a  cause  or  an  effect  of  crack  behavior,  the  measurements  were 
useful  in  studying  changes  at  the  crack  tip  during  holding  at 
zero  load. 

Conclusions 

1.  As  hold  time  at  zero  load  increased  the  number  of  delay 
cycles  decreased  for  7075-T651,  2219-T851  and  overaged 
7075  and  2219  alloys  in  humid  air.  This  was  also  true 
for  7075-T651  alloy  in  vacuum.  This  may  be  due  to  a 
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time  dependent  strain  relaxation  of  the  residual  strain 
created  by  the  overload. 

2.  The  relaxation  effects  were  greater  in  humid  air  than  in 
vacuum. 

3.  For  the  same  load  spectra,  the  overaged  alloys  exhibited 
shorter  retardation  times  than  high  strength  alloys. 

4.  The  delay  periods  for  7075-T651  in  vacuum  were  signifi¬ 
cantly  higher  than  in  humid  air. 

5.  Changes  in  crack  closure  during  retardation  indicate  that 
surface  hammering  may  be  occurring  during  retardation. 
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A  MODEL  FOR  FATIGUE  CRACK  CLOSURE  BASED  OH 
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Introduction 

The  crack  closure  aodel,  first  developed  by  Elber  (1),  attempts  to  explain  spectrum- 
loading  fatigue  crack  growth  in  terns  of  the  effective  stress  intensity  range, OK  ... 
According  to  this  aodel,  as  a  fatigue  crack  grows,  it  leaves  behind  a  wake  of  res Huai  strain 
(Figure  1)  which  aay  cause  the  crack  surfaces  to  touch  during  the  lower  part  of  a  tension- 
tension  loading  cycle.  This  crack  closure  has  been  observed  experimentally  (1,2, 3,4),  It 
is  thought  that  once  the  eurfaces  touch  there  is  no  longer  a  singularity  in  the  stresses  st 
the  crack  tip,  i.e.,  the  stress  intensity  is  aero.  Thus  the  effective  stress  intensity 
range,  which  determines  the  crack  growth  rate,  should  be  calculated  fros  the  effective  load 
range,  that  is,  the  range  of  loads  for  which  the  crack  is  open,  so  that  AK 
where  P  is  the  load  st  which  the  crack  opens  and  C  Is  the  geoaetrleal 
stress  Intensity  to  load.  C  la  s  function  of  the  craok  length. 

Coaputer  Modal  for  Crack  Closure 

In  this  section  we  develop  a  coaputer  aodel  for  crack  closure  during  constant  aaplitude 
loading.  The  coaputer  aodel  will  be  used  to  predict  the  opening  load  and  the  extent  of 
closure  at  lower  loads.  The  comparison  of  these  predictions  with  experiment  is  given  in 
the  next  section. 

Aasuae  that  a  fatigue  crack  growth  experlaent  Is  carried  out  on  a  cospact  tension 
specimen  with  loading  controlled  so  as  to  aalntaln  constant  K  .  and  K.  As  the  crack 
grows,  it  will  leave  behind  an  approxlaately  uniform  wake  of  residual  strain,  that  is,  tha 
■ate rial  near  the  crack  surfaces  will  effectively  becoae  thicker,  because  of  the  surface 
roughness  and  voluae  change  associated  with  plastic  defoliation  by  high  tensile  stresses 
near  the  crack  tip  before  final  separation.  The  Increase  in  the  thickness  can  be  modeled  aa 
due  to  a  uniform  coating  of  thickness  d  of  new  material  on  the  crack  surfaces  (Figure  1), 
Extending  in  front  of  the  crack  tip  (B)  for  a  distance  equal  to  the  plastic  sone  sice  r -BC , 
there  will  also  be  "extra  material",  decreasing  in  thickness  from  2d  st  the  tip  to  cero” 
outside  the  plastic  cons.  If  this  "extra  material"  were  cleared  away  the  new  crack  surfaces 
(we  will  call  than  the  "aodel  crack"  surfaces)  would  just  touch  without  exerting  any  forces 
on  each  other  when  the  external  load,  P,  is  resoved. 

Thus  our  aodel  consists  of  a  ssapls  with  a  "aodel  crack",  longer  than  the  true  crack 
by  a  distance  r  ,  with  some  incompressible  external  material  Inserted  between  the  crack 
surfaces  propping  them  open.  The  true  crack  say  actually  be  closed  a  distance  AB  behind  the 
true  crack  tip,  with  AB  depending  on  the  external  load. 

For  a  given  load,  P,  we  use  a  flnlte-eleaent  elastic  stress  analysis  program  written  by 
one  of  the  authors  ( 5)  to  find  the  stress  state  of  the  saaple.  One  of  the  aost  important 
features  of  this  program,  which  will  be  available  through  the  slgor lthas  distribution  service 


i'ATIfiUf:  CRACK  cinsuiu; 


Vol .  11,  No.  b 


of  (5)*  i*  ability  to  efficiently  handle  crack  problems.  If  we  knew  the  closure  distance 
beforehand,  we  could  give  the  following  boundary  conditions  along  the  aodel  crack  surfaces i 
Free  the  point  of  closure  (a)  on  out  away  fros  the  crack  tip,  we  know  that  the  surfaces  are 
free  surfaces.  Fros  the  point  of  closure  on  In  towards  the  crack  tip,  we  know  the 
displacements ,  naaely,  from  A  to  B  the  surfaces  are  displaced  vertically  by  a  distance  d, 
and  froa  B  to  C  the  displacements  decrease  from  d  to  0.  In  our  aodel  they  decrease  as  a 
cubic  polynomial  with  sero  slope  at  B  and  C,  i.e. ,  as  d(3s2-2s3)  where  e*(C-x)/(C-B). 

Unfortunately  we  cannot  determine  the  closure  distance  AB  before  solving  the  problem, 
so  we  proceed  by  an  iterative  nethod.  If  our  guess  for  AB  is  short,  this  means  we  are 
allowing  part  of  the  crack  surfaces  to  be  free  which  would  in  reality  be  touching,  so  our 
resulting  calculations  will  show  that  the  true  crack  surfaces  overlap  behind  the  assumed 
point  of  closure.  If  our  guess  is  long,  this  means  we  are  glueing  part  of  the  surfaces 
together  that  would  in  reality  not  be  touching,  so  that  the  true  crack  will  be  open  behind 
the  assumed  point  of  closure.  The  assumed  point  of  closure  will  be  In  compression  or  tension 
depending  on  whether  the  guess  is  short  or  long,  respectively,  but  the  stresses  are  not 
calculated  ae  accurately  as  the  displacements,  so  we  do  not  use  them  to  guide  the  iteration. 

A  series  of  computations  were  made  with  the  program,  for  a  compact  tension  sample  with 
W-0.108*,  H-0.053m,  thlckness-0. 025a,  EWOxloSw/m2  and  v-0.33  with  plane  etrmin  conditions. 
The  true  crack  length  was  taken  as  a-0.056a  with  a  plastic  sone  sire  of  r  -0.0014a  and  s 
value  for  d  of  25  u-n  was  used.  The  loads  corresponding  to  three  closure  distances  are  found 
by  iteration  and  tabulated  below. 


TABLE  1 

Loads  for  Three  Closure  Distances  (d-BS^s) 

Closure  (AB)  Load  (P) 

0.0111m  -1*5  kg 

0.0056m  320  kg 

0  (Just  open)  1840  kg 


The  linearity  of  the  problem  implies  that  the  values  of  d  and  P  may  be  multiplied  by  a 
common  factor  without  changing  the  closure  distance.  Thus  for  this  problem,  with  the  given 
crack  length,  our  aodel  predicts  an  opening  load  of  ?.23xl07»d  kg,  if  d  is  in  meters.  Also, 
by  Interpolation,  we  see  that  the  closure  distance  ehould  be  about  0.010a  for  rero  applied 
load,  regardless  of  the  value  of  d. 

The  vertical  stress  distribution  along  the  closure  surfaces  and  Into  the  plastic  sons 
Is  shown  for  P-  -45  kg  In  Figure  2. 

Comparison  of  Coaputer  Simulation  with  Experimental  Results 

To  check  the  crack  closure  predictions  given  In  the  previous  section,  two  A1-2219-T851 
compact  tension  specimens  of  the  same  dimensions  and  elastic  properties  as  In  the  coaputer 
simulation  ware  tested  on  an  KTS  closed-loop  electrohydraullc  testing  machine,  after 
overaging  for  20  hours  at  J00°C  In  air.  Tensile  specimens  given  the  esse  heat  treatment 
wera  found  to  have  yield  strengths  of  about  200  MN/m^.  One  purpose  of  the  lowering  of  the 
yield  strength  was  to  give  a  larger  and  more  easily  neasureable  value  of  d,  the  residual 
strain  parameter. 

In  each  of  the  two  samples  tested,  a  fatigue  crack  was  grown  to  a-0,056a  with  loading 
controlled  so  that  K  .  -0  and  were  held  approximately  constant.  An  extensiometer  was 
placed  so  as  to  straddle  the  craw  (at  a-0 , 044« ,  point  C  of  Figure  3)  as  it  was  grown 
through.  The  Increase  In  distance  between  the  legs  of  the  extensloaeter  as  the  crack  passed 
between  them  (distance  always  measured  with  the  crack  closed)  Is  taken  to  be  the  residual 
strain,  2d  (Figure  4),  since  the  legs  are  far  enough  apart  to  straddle  all  the  plastic 
strain  and  crack  front  roughness.  To  determine  the  opening  load,  P00#n»  and  the  loede  P_ 
and  P„  at  which  the  crack  closes  to  points  B  and  C  of  Figure  3  (representing  closure  of 
0.005oa  and  0.0111m  respectively),  ultrasonic  transducers  ware  placed  in  succession  at 
positions  1,  2  and  3,  The  two  piezoelectric  transducers  transmit  and  receive  a  2.25  mega- 
herts  longitudinal  ultrasonic  wave  which  Is  converted  to  s  DC  signal.  This  signal  Is 
attenuated  as  the  crack  grows  through  the  path  of  the  wave  (see  reference  4),  so  that  in 
position  1,  the  signal  increases  when  the  crack  first  begins  to  close  anl  in  positions  2 
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or  3  the  signal  la  zero  until  the  crack  closes  to  point  B  or  C  respectively.  Figure  $  shows 
the  signal  as  a  function  of  load  for  each  position  of  the  transducer.  Both  Figures  4  and  5 
are  copied  directly  fro®  the  X-Y  recorder  output. 

Loading  of  the  first  sample  was  a  2-hertz  sine  wave  with  K  -26.4  MPa /S'.  The  plastic 
cone  alee  is  estiaated  as  r  »(K  ^9'y)2/4 JZrr  -0.0010a  (6),  The  value  of  d  measured  by  the 

extensiometer  was  8,1^ a,  wnichDa&ee  it  possible  to  calculate  the  eoaputer  predictions  for 

Popen*  PB  and  PC  ** 

The  second  sample  was  loaded  with  a  5-hertz  sine  wave  with  K-13.2  MPa/a,  and  d  was 
measured  as  2.5/< The  experimental  results  for  both  samples  arecoapared  with  the 
corresponding  computer  predictions  below,  Ignoring  the  fact  that  the  plastic  zone  size  for 
the  second  saaple  is  four  times  smaller  than  for  the  first  sample. 


TABLE  2 

Experimental  and  Calculated  Closure  Loads 


First  Sample 

P 

open 

PB 

PC 

Second  Sample 

P 

open 


Experimental  (+50  kg) 

556  kg  " 

100  kg 
0 

Experimental  (+20  kg) 
150  kg 

50  kg 
-100  kg 


Calculated  (d-8.lz.rn) 
590  kg 

100  kg 

-15  kg 

Calculated  (d-2«5/*a) 
180  kg 

30  kg 
-5  kg 


Discussion 

The  computer  predictions  of  the  opening  load  and  extent  of  closure  at  lower  loads  for 
s  crack  In  a  compact  tension  specimen  showed  close  agreement  with  the  experimental 
measurement. 

The  elastic  analysis  shows  that  the  opening  load  Increases  as  d,  a  measure  of  the 
residual  strain  and  crack  surface  roughness,  Increases.  This  Is  consistent  with  the  crack 
closure  model  explanation  of  the  well-known  effects  of  overloads  (3, 4, 7, 8,9)  and  underloads 
(7,10).  The  effect  of  an  overload  Is  to  generate  a  large  amount  of  residual  strain  (3,4) 
which  causes  the  crack  to  close  st  a  higher  load  when  baseline  loading  la  resumed  and  thus 
dK  Is  decreased  and  crack  growth  is  retarded,  A  subsequent  underload  has  the  opposite 
effect,  by  decreasing  the  residual  strain  through  compression,  and  thus  increasing  AK 
and  the  crack  growth  rate. 
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Introduction 

During  fatigue  crack  growth,  plastic  deformation  is  left 
in  the  wake  of  the  crack.  As  a  result  of  this  residual  defor¬ 
mation,  the  crack  can  physically  close  during  the  lower  part 
of  the  tensile  crack  unloading.  This  closure  phenomenon  has 
been  observed  and  studied  extensively^1-10^.  Various  measure¬ 
ment  techniques  have  been  employed  to  measure  crack  closure, 
and  many  interesting  observations  are  presented  on  the  basis 

of  those  measurements.  With  the  help  of  potential  drop  measure- 

1’51 

ments  Lindley  and  Richards  concluded  that  crack  closure  is 

only  a  surface  or  edge  effect.  Potential  drop  measurements 

have  been  used  by  other  researchers^'10^.  Irving^  did 

experiments  in  vacuum  as  well  as  humid  air  and  observed  crack 

closure  only  in  vacuum.  Bachman  and  Munz^10^  however,  found 

that  closure  is  present  in  vacuum  as  well  as  at  760  Torr  pressure 

and  that  the  crack  opening  stress  intensity  (K  )  was  indepen- 

[12] 

dent  of  environment.  Avery  important  conclusion  of  their 
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study  was  that  in  humid  air  environment,  an  oxide  layer 

develops  on  the  crack  surfaces,  and  this  effect  might  give  some 

erroneous  results  in  terms  of  the  potential  drop  measurements. 

They  point  out  that  earlier  results  based  on  the  potential  drop 

measurements  can  be  reinterpreted. 

r  7 1 

Mahulikar  et  a_l  have  discussed  edge  effects  on  crack 
closure.  Because  of  the  plane  stress  region  at  the  edges  of 
the  compact  tension  specimen, surface  measurements  are  not  repre¬ 
sentative  of  the  bulk  behavior.  Ultrasonic  bulk  measurements 

have  demonstrated  that  crack  closure  is  a  plane  strain  as  well 

f  7  31 

as  a  plane  stress  effect  ' 

Studies  of  the  residual  displacements  left  in  the  wake  of 
the  crack,  can  give  some  important  information  about  crack  clo¬ 
sure.  The  residual  displacement  is  considered  to  be  made  up 
of  the  nonuniform  plastic  strain  on  the  surface,  the  small 
change  in  volume  due  to  the  deformation  process  and  the  oxide 
formed  on  the  nonuniformly  strained  fracture  surface.  It  was 
demonstrated^11-*  that  the  residual  displacements  can  be  measured 
by  a  clip  on  gage  mounted  on  the  side  of  the  specimen.  Sewell 
and  Marcus ^11^  measured  the  residual  displacements,  (d) ,  for 
different  specimens  in  humid  air  and  proposed  a  finite  element 
model  to  predict  crack  opening  load  (P  )  using  the  ’d’  measure¬ 
ment. 

In  order  to  clearly  define  the  relationship  of  d  with 
PQp,  it  was  necessary  to  perform  experiments  in  vacuum  as  well 
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as  other  environments.  Contributions  of  plane  stress  region 
and  oxidation  to  the  residual  displacement  must  be  studied. 

This  paper  presents  the  results  for  such  experiments  on  alumi¬ 
num  alloys.  An  effort  has  been  made  in  this  paper  to  relate 
ultrasonic  bulk  measurements  of  crack  closure  to  the  surface 
measurements  of  the  residual  displacement.  The  results  clearly 
demonstrated  that  surface  measurements  have  limitations  and  use 
of  such  surface  measurements  for  predicting  certain  aspects  of 
fatigue  crack  growth  can  be  misleading. 

Experimental  Procedure 

Experiments  were  performed  on  the  7075-T651  and  2219-T851 

aluminum  alloys.  Compact  tension  specimens  with  dimensions 

B  =  2.54  cm  (thickness),  H  =  5.3  cm  (halfheight)  and  VJ  =  10.8 

cm  (width)  were  used.  Four  types  of  experiments  were  performed. 

In  a  first  set  of  experiments  only  one  type  of  environment 

(e.g.  vacuum)  was  used.  The  crack  was  initiated  and  grown  in 

that  environment  at  constant  AK.  Closure  curves  were  obtained 

f  3] 

by  ultrasonic  techniques  described  elsewhere  .  Residual  dis¬ 
placements  were  measured  by  a  clip  on  gage  mounted  on  the  side 
of  the  specimen.  Since  the  legs  of  the  clip  on  gage  are  approxi¬ 
mately  7  mm  in  length,  it  was  necessary  to  take  two  readings  to 
determine  the  residual  displacement  d  .  The  initial  reading 
(Fig.  2)  was  taken  just  before  crack  entered  the  gage  and  the 
final  reading  was  taken  just  after  the  crack  left  the  gage. 

The  technique  to  obtain  d  value  is  described  in  [11] .  The 
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geometry  of  the  combined  measurement  is  shown  in  Fig  1.  The 
gages  were  mounted  at  a/w  =  0.48. 

In  another  set  of  experiments  the  crack  was  first  grown 
in  vacuum  and  the  'd'  measurement  curve  was  obtained  by  the 
gage  mounted  at  a/w  =  0.48.  The  cycling  was  then  stopped  and 
the  specimen  held  at  the  mean  load  (P  ) .  D-0  water  environ- 

ment  (=100%  RH)  was  then  introduced.  As  soon  as  the  D2O  water 
was  introduced  the  residual  displacement,  d  ,  was  obtained. 
Cycling  was  then  continued  at  constant  AK,  and  after  1000  cycles 
the  d  measurement  was  repeated. 

In  the  third  set  of  experiments,  the  crack  was  grown  in 
vacuum,  but  before  introducing  D20  water,  dry  nitrogen  gas  was 
introduced  and  the  d  curve  was  obtained.  D20  water  was  then 
introduced  and  d  measurement  repeated.  The  fourth  type  of 
experiment  on  the  7075  alloy  consisted  of  introducing  dry  C>2 
gas  after  the  introduction  of  N2  gas. 

For  ultrasonic  measurements,  two  2.25  MHz  longitudinal 
transducers  of  12.7  cm  dia.  were  used.  All  the  experiments 
were  done  in  an  all  metal,  ion  pumped,  environment  chamber 
mounted  on  an  MTS  electrohydraulic  testing  machine.  The  fre¬ 
quency  of  cycling  during  fatigue  crack  growth  was  6  Hz.  The 
ultrasonic  and  d  measurements  were  made  at. 06  Hz.  Loading 
parameters  for  each  experiment  are  presented  with  the  results. 

Results 

Fig.  2  and  3  show  typical  ultrasonic  transmission  curves 
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as  a  function  of  applied  load  for  closure  measurement  and  d 
measurements  respectively.  In  Fig.  1,  represents  the  load 
at  which  the  curve  deviates  from  linearity,  while  P2  represents 
the  closure  load  obtained  from  the  tangent  method  In  Fig.  2 

the  initial  d  displacement  curve  represents  the  reading  taken 
before  the  crack  enters  the  gage  and  the  final  curve  represents 
the  reading  taken  just  after  the  crack  leaves  the  gage.  The 
x-axis  shift  represents  2  d  units  since  d  is  defined  as  the 
residual  displacement  on  one  fracture  surface. 

Fig.  4  shows  the  effect  of  environments  on  the  closure 
load  (P2)  and  on  the  d  values  observed  for  7075-T651  alloys. 
Fig.  5  shows  the  same  effect  for  the  2219-T851  specimens. 

Fig.  6  plots  the  crack  growth  rates  against  the  effective  stress 
intensity  AKe££  as  calculated  from  P2  for  both  the  alloys. 

The  d  curves  presented  in  Fig.  7  are  for  the  four  types 
of  experiments  defined  in  the  experimental  procedure  section. 

Fig.  7  (a)  is  representative  of  the  first  type  of  experiment; 

Fig.  7  (b)  represents  the  third  type  of  experiment  and  Fig.  7 
(c)  represents  the  fourth  type  of  experiment. 

Discussion 

The  following  observations  can  be  made  from  Figs.  3  and  4. 

a)  The  closure  load  appears  to  decrease  with  increase  in 
humidity.  It  also  decreases  as  the  d  value  increases. 

b)  The  value  of  d  increases  with  the  introduction  of  humidity. 

c)  The  values  of  d  for  the  higher  strength  7075-T651  alloy 
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are  higher  than  the  2219-T851  alloy;  and  so  are  the  closure 
loads . 

Thus,  it  was  observed  that  the  closure  load  bears  an 
inverse  relationship  with  the  d  values.  This  unexpected 
result  can  be  explained  by  the  fact  that  d  measurement  is  a 
surface  measurement  and  ultrasonic  measurement  of  closure  is 
a  bulk  measurement.  It  should  be  noted  that  the  ultrasonic 
transducers  of  1.27  cm  dia.  were  mounted  on  a  specimen  2.54 
cm  thick  at  the  middle  portion.  Thus  they  were  sampling  a 
columnar  region  of  approximately  1.27  cm  dia.  in  the  center 
of  the  specimen.  The  plane  stress  edge  region  of  the  speci¬ 
men  was  not  being  sampled.  On  the  edges,  the  residual  dis¬ 
placement  was  higher  than  the  region  sampled  by  the  trans¬ 
ducers.  Thus  even  though  closure  occurred  first  on  the  edges 
the  ultrasonics  would  not  detect  it  until  the  inner  portion 
and  largest  portion  of  the  crack  front  was  compressed  suffi¬ 
ciently  to  touch  the  surfaces.  On  the  other  hand,  the  d 
measurements  shown  are  a  surface  measurement  representative 
of  the  plane  stress  plastic  deformation  at  a  distance  behind 
the  crack  tip.  In  humid  air  environment  a  faster  buildup  of 
oxide  layer  on  the  surfaces  of  the  crack  occurs.  This  effect 
can  be  expected  to  be  more  severe  on  the  plane  stress  surface 
region  of  the  specimen.  The  'extra'  material  consisting  of 
the  oxide  layer  as  well  as  mechanical  strengthening  of  asperi 
ties  by  the  oxide  on  the  edges  of  the  crack  front  results  in 
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increased  d  value  in  humid  air  beyond  that  of  just  the  oxide 
thickness.  Therefore  during  the  reverse  loading  more  material 
is  available  for  contact  on  the  edges  of  the  specimen  than  in 
the  middle  of  the  specimen.  Therefore  the  crack  front  in  the 
middle  of  the  sepcimen  can  only  touch  after  the  load  is  reduced 
further  and  a  lower  closure  load  is  observed  on  the  ultrasonic 
signal. 

For  7075  alloy  the  d  value  in  dry  0^  gas  is  lower  than 
for  humid  air  as  would  be  expected  due  to  the  slower  formation 
of  the  oxide  layer  in  dry  C>2  gas  than  in  humid  air .  The  crack 
growth  rate  is  seen  to  increase  with  the  humidity  (Fig.  6)  as 
does  the  AKeff  calculated  using  P,,. 

An  interesting  observation  was  made  during  the  second  set 
of  experiments  (Fig.  7a) .  The  crack  was  grown  through  the  clip 
on  gage  in  vacuum  and  the  d  curve  was  obtained  as  usual.  The 
cycling  was  then  stopped  and  the  specimen  was  held  at  the  mean 
load.  D2O  water  vapor  to  100%  RH  was  then  introduced  in  the 
environment  chamber.  As  soon  as  the  water  vapor  was  introduced, 
the  gage  recorded  a  shift  (Fig.  7a)  in  the  displacement  reading. 
This  shift  was  of  the  order  of  3  to  10  ym  for  the  7075  alloys. 

The  shift  increased  slowly  with  time.  After  about  an  hour  it 
was  about  12  ym.  After  this  cycling  was  continued  for  1000 
cycles  and  a  d  curve  was  recorded  with  a  smaller  further  incre¬ 
ment  to  the  value  of  d.  This  shift  in  d  value  after  the  intro¬ 
duction  of  the  D20  water  vapor  is  most  likely  due  to  the  formation 
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of  oxide  layer  which  gives  mechanical  strengthening  of  the 
asperities  on  the  edge  of  the  specimen.  However,  the  magni¬ 
tude  of  this  shift  is  certainly  not  representative  of  oxide 
layer  thickness  on  aluminum  surfaces  which  leads  one  to  believe 
that  some  additional  mechanisms  associated  with  the  oxide  for¬ 
mation  must  be  operating.  The  strengthening  of  the  asperities 
is  a  possible  mechanism. 

In  order  to  establish  that  the  change  in  pressure  in  the 
system  was  not  the  cause  of  the  d  shift,  the  third  set  of 
experiments  were  carried  out.  After  the  crack  was  grown  in 
vacuum  but  before  introducing  D2O  water,  dry  N2  gas  was  intro¬ 
duced  (Fig.  7b).  This  showed  no  change  or  shift  in  'd'  reading. 
This  eliminated  the  possibility  of  a  strictly  mechanical  effect 
from  the  pressure  change  to  explain  the  change  in  d.  D20  water 
was  then  introduced,  which  again  showed  a  significant  shift  in 
the  d  value  that  increased  even  further  on  subsequent  cycling. 

The  experiment  was  then  repeated  with  dry  02  gas  intro¬ 
duced  after  the  vacuum  and  N2>  The  shift  in  d  was  very  small 
'1  pm  and  further  cycling  in  the  dry  02  showed  only  a  small 
increase  in  d  value,  Fig.  7c. 

Thus  it  was  concluded  that  oxidation  played  a  dominant 
role  in  this  shift.  Oxidation  in  a  dry  02  gas  environment  is 
a  very  slow  process  when  compared  to  oxidation  in  a  humid 
environment.  This  fast  oxidation  in  a  humid  environment  must 
attack  the  plane  stress  surface  portion  of  the  crack 
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preferentially.  This  results  in  changes  in  the  residual  dis¬ 
placement  probably  due  to  mechanical  strengthening  by  the  oxides 
of  asperities  that  are  more  dominant  on  the  edge  of  the  speci¬ 
men.  These  effects  yield  the  shift  in  d  readings,  while  the 
edge  effect  holding  the  crack  open  results  in  a  lower  closure 
load  in  the  center. 

It  was  thus  apparent  that  a  lot  of  changes  on  the  edges 
or  surfaces  of  the  crack  front  take  place  when  one  goes  from 
vacuum  to  humid  air.  Surface  readings  can  be  entirely  different 
than  the  bulk  readings  because  of  effects  of  oxidation,  such 
as  the  strengthening  of  asperities.  Surface  clip  on  gage 
measurements,  thus,  cannot  be  considered  totally  representative 
of  the  bulk  closure  behavior.  Thus  the  use  of  d  surface  measure¬ 
ment  for  load  predictions  in  a  varying  environment  must  be  care¬ 
fully  evaluated.  In  fact  using  the  data  presented  here,  the  d 
measurements  combined  with  a  finite  element  analysis  would  pre¬ 
dict  a  higher  crack  closure  load,  a  lower  and  reduced 

growth  which  is  contradictory  to  experimental  observations. 

Conclusions 

1.  In  7075-T651  and  2219-T851  the  closure  load  decreased  with 
the  introduction  of  humid  environments  for  both  the  alloys, 

while  residual  displacement  d  values  increased, 

2.  The  decrease  in  the  closure  load  with  humidity  is  attributed 
to  buildup  of  extra  material  at  the  edges  of  the  specimen 
which  the  acoustic  transducers  do  not  sample. 
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3.  The  increase  in  d  value  when  fatigued  in  humid  environ¬ 
ment  compared  to  vacuum  is  attributed  to  the  buildup  of 
oxides  on  the  edges  of  the  fracture  surface. 

4.  The  d  value  was  observed  to  increase  when  D20  water 
environment  was  introduced  after  the  crack  was  grown  in 
vacuum.  Dry  02  gas  introduced  after  crack  growth  in 
vacuum  showed  only  a  very  small  shift  in  d  value. 

5.  The  shift  in  d  value  resulting  from  the  introduction  of 
D20  water  vapor  after  the  crack  was  grown  in  vacuum  was 
found  to  increase  with  time  and  additional  cycling. 

6.  It  is  suggested  that  in  addition  to  oxide  formation  itself, 
preferential  strengthening  of  asperities  at  the  edge  of 
the  specimen  may  explain  the  increase  in  d  in  humid  air 
environments . 

7.  Surface  measurements  are  not  always  representative  of  the 
bulk  behavior  and  use  of  such  measurements  in  predictive 
analysis  can  be  misleading. 
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Int  roduct ion 


The  crack  closure  effect  suggested  by  fiber  (1)  has  been  studied  quite 
extensively.  The  fact  that  the  crack  closes  in  the  lower  part  of  load  cycle 
has  been  observed  and  measured  universally  but  the  extent  to  which  the  closure 
affects  retardation  phenomenon  is  not  yet  known  clearly.  There  have  been 
speculations  that  the  crack  closure  effect  is  only  a  surface  or  edge  effect. 
This  was  supported  by  experiments  of  Lindley  and  Richards  (2).  They  used 
potential  drop  measurements  as  well  as  metallographic  techniques  to  study 
closure  effects  in  the  steels.  They  concluded  that  crack  closure  is  only  a 
surface  or  edge  effect.  This  possibility  is  also  discussed  by  other  research¬ 
ers  (3,4).  It  was  suggested  that  the  crack  closes  only  at  the  surface  or  edges 
of  the  specimen  where  plane  stress  conditions  exist,  while  it  may  not  close  In 
the  central  part  where  plane  strain  conditions  exist.  This  paper  reports  on 
experiments  aimed  at  clarifying  this  point. 

Experimental  Procedures 

Several  experiments  were  performed  to  study  this  edge  effect. 
Experiments  were  done  on  707S-T651,  2219-T8S1  and  overaged  2219  aluminium 
alloys.  The  mechanical  properties  of  these  materials  are  listed  in  Table  1. 
Compact  tension  specimen  with  dimensions  B*2.54  cm  thick,  H*5.3  cm  halfheight 
and  W-10.8  cm  long  were  used. 

Two  types  of  experiments  were  performed  on  the  7075-T651,  2.54  cm  thick 
specimens.  In  the  first  set,  the  crack  was  grown  up  to  4.31  cm  at  a  constant 
AK  (16.5  MPa/m)  and  the  load  ratio,  R-0.17,  on  a  closed  loop  hydraulic  testing 
machine.  Closure  measurement  was  made  by  two  1.27  cm  diameter  ultrasonic 
transducers  mounted  directly  above  and  below  the  crack  tip  and  centered  on  the 
width  of  the  specimen.  The  ultrasonic  technique  of  load  and  crack  length 
measurement  is  described  elsewhere  (6).  In  our  experiments,  the  ultrasonic 
wavelength  was  significantly  longer  than  the  crack  tip  radius,  hence  the  pos¬ 
sibility  of  closure  curve  changes  being  only  diffraction  or  scattering  effects 
was  considered  unlikely.  The  experiment  was  then  stopped  and  two  sides  of  the 
specimen  were  machined  down  to  a  thickness  of  1.75  cm.  Again  the  closure 
readings  were  taken  at  the  same  point.  The  cycling  was  resumed  at  the  same  AK 
value  and  the  crack  was  grown  to  about  4.8  cm  and  closure  reading  were  taken 
directly  above  the  crack  tip.  The  same  procedure  was  repeated  by  machining 
down  the  specimen  to  a  thickness  of  1.25  cm,  taking  closure  measurement  at  that 
point,  and  again  after  growing  the  crack  to  a  length  of  6.3  cm.  The  results 
are  presented  in  Fig.  1  where  the  applied  stress  intensity  is  plotted  against 
t he  received  acoustic  intensity. 
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In  a  second  set  of  experiments  a  specimen  used  for  relaxation  experi¬ 
ments  (5)  was  used  for  edge  effect  experiments.  A  fatigue  crack  was  grown  in 
a  707S-T6S1  specimen  to  a  length  of  4.84  cm  at  a  AK  of  6.6  MPa/m,  R*0.33. 

It  was  allowed  to  rest  at  a  zero  load  for  a  period  of  77  hours  for  relaxation 
measurements.  Closure  measurements  were  made  at  the  end  of  the  zero  load  hold 
time.  Then  the  specimen  was  machined  down  to  a  thickness  of  1.3  cm,  thus 
insuring  that  entire  plane  stress  region  was  removed.  Closure  measurements 
were  made  at  this  point.  Cycling  was  resumed  and  the  crack  was  grown  to 
around  S.4  cm  length  with  continual  measurement  of  closure. 

Similar  tests  were  performed  on  2219-T851  and  overaged  2219  alloys. 

In  these  tests  the  crack  was  grown  at  a  frequency  of  2  llz ,  AK*26.4  UPa/m,  and 
R-0  to  a  length  of  S.6  cm.  At  this  point  cycling  was  stopped  and  closure 
measurements  were  made.  Both  sides  of  specimen  were  then  machined  off  to  a 
thickness  of  1.9  cm  and  closure  was  measured  with  the  transducers  at  same 
position  over  the  crack  tip.  The  specimen  was  then  machined  off  from  both 
sides  to  a  thickness  of  1.3  cm  and  once  more  closure  measurements  were  made. 

It  should  be  noted  that  the  original  specimen  thickness  was  2.54  cm, 
twice  the  diameter  of  the  1.27  cm  ultrasonic  transducers.  The  ultrasonic  signal 
outside  this  1.27  cm  columnar  region  does  not  contribute  significantly  to  the 
amplitude  of  the  output  because  of  the  method  used  in  converting  the  signal  in 
the  peak  to  peak  detector.  At  a  specimen  thickness  of  2.54  cm  the  signal 
transmitted  through  the  plane  stress  region  does  not  contribute  significantly 
to  the  output  signal  intensity  since  the  transducers  are  centered  across  the 
thickness  of  the  specimen. 

Results  and  Discussion 

Closure  curves  are  presented  in  Fig.  1.  Two  stress  intensities  Kj  and 
Kj  were  obtained  from  the  curves.  Kj  is  the  intensity  at  which  the  closure 
curve  deviates  from  linearity  while  K?  is  the  closure  intensity  determined  by 
the  tangent  method  described  by  Buck  et"al.  (6).  Tables  2,  3,  4  and  5  list 
all  the  values. 

From  the  ultrasonic  measurements  for  both  the  alloys  and  all  cases, 
significant  crack  closure  was  still  observed  after  machining  off  the  plane 
stress  region  as  is  evident  from  Fig.  1  and  Tables  2,  3,  4  and  5. 

For  7075  alloys,  it  was  observed  that  machining  off  the  sides  did  not 
change  the  closure  stress  intensities  significantly.  However,  as  the  crack 
was  grown  further,  the  values  of  and  Kj  increased.  On  the  other  hand  for 
2219  alloys  a  significant  increase  in  Kj  and  was  observed  after  machining 
off  the  siJcs. 

These  observations  can  be  explained  from  the  fact  that  7075  is  a  high- 
strength,  low  Kjc  alloy.  Hence  the  plastic  zone  sizes  encountered  are 
relatively  smaller.  Thus  within  the  limits  of  ultrasonic  measurements  Kj ,  K2 
did  not  change  significantly.  When  the  crack  was  grown  further,  new  plane 
stress  region  developed  indicating  slight  rise  in  Ki  and  K2  at  a  thickness  of 
1.75  cm  and  a  significant  rise  at  a  thickness  of  1.25  cm.  The  significant 
rise  can  be  understood  by  the  fact  that  the  thickness  is  very  small,  slightly 
greater  than  allowable  thickness  for  plane  strain  testing,  and  transducers  are 
now  sampling  the  entire  thickness  of  the  specimen.  Higher  AK  values  were  not 
used  due  to  the  limitations  of  the  fracture  toughness. 

For  2219  alloys,  machining  off  the  sides  of  t he  specimens  instanta¬ 
neously  increased  the  values  of  Kj  and  K7 .  ZZ19  has  a  lower  Oy  and  higher  Kjc 
than  7075.  The  higher  made  it  possible  to  use  higher  AK  values  for 
the  experiment.  This  combined  with  the  lower  value  of  Oy  makes  the  plastic 
zone  sizes  significantly  larger.  When  the  specimen  was  at  the  original 
thickness  this  large  plastic  strain  in  the  plane  stress  region  allowed  the 
specimen  edges  to  close  before  the  plane  strain  region.  This  served  to  prop 
open  the  interior  and  keep  it  from  closing  until  the  load  was  reduced  enough 
to  sufficiently  compress  the  plane  stress  region  so  the  interior  closed.  The 
result  was  a  lower  closure  stress  intensity  detected  by  ultrasonics.  Once  the 
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plane  stress  region  was  machined  off  the  remaining  plane  strain  region  would 
close  at  a  higher  value  of  K. 

A  brief  mention  should  be  made  of  the  residual  stresses.  Residual 
stresses  are  introduced  on  the  surface  of  the  specimens  because  of  heat  treat¬ 
ments  given  to  them  and  also  because  of  machining.  The  extent  to  which  these 
stresses  can  affect  the  closure  behaviour  is  not  yet  well  understood,  but  are 
not  expected  to  be  significant  in  these  experiments. 

Conclusions 

1.  Significant  crack  closure  was  ohserved  for  ultrasonic  measurements 
made  after  machining  off  the  plane  stress  region  of  the  fatigue  crack  in  7075 
and  2219  aluminium  alloys. 

2.  For  high  strength  low  Xjc  7075  alloys  closure  stress  intensities 
did  not  change  significantly  after  machining  off  the  sides  of  the  compact 
tension  specimen,  however  Kj  and  K,  increased  after  the  crack  was  grown  further. 

3.  For  relatively  low  strength  high  Kic  2219  alloys  the  value  of  Ki  and 
K2  increased  after  machining  off  the  sides.  This  is  attributed  to  relatively 
large  plastic  zone  sizes  in  2219. 

4.  Crack  closure  is  not  only  associated  with  the  edge  of  the  specimen, 
but  with  the  total  crack  front  in  compact  tension  specimens. 
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TABLE  1 


Yield  Strengths  and  Kj  Values  of  Specimen  Materials 


Material 

UyMPa 

K](.MPa/m 

7075-T651 

505 

25 

2 219-T851 

350 

35 

2219-overaged 

200 

36 

TABLE  2 

Material:  7075-T651;  AK 

=  16.5  MPa/in 

Thickness  Crack  Length  K, 

(MPa/m) 

K2  (MPa/ii) 

B  (cm) 

a  (cm) 

2.54 

4.3 

8.3 

5.005 

1.75 

4.3 

7.3 

4 .334 

1.75 

4.8 

8.2 

5.225 

1.25 

4.8 

8.35 

4.455 

1.25 

6.3 

TABLE  3 

16.8 

8.36 

Material:  7075  T6S1;  AK 

=6.6  MPa/m 

B  (cm) 

a  (cm) 

K1 

(MPa/in) 

K2  (MPa/in) 

2.5 

4.8 

2.43 

1.73 

1.3 

4.8 

3.31 

1.56 

1.3 

5.4 

TABLE  4 

3.34 

1.95 

Material:  2129-T8S1;  AK 

=  26.4  MPa /in 

B  (cm) 

a  (cm) 

K1 

(MPa/iff) 

K2  (MPa/in) 

2.5 

5.6 

7.6 

4.6 

1.9 

5.6 

15.6 

5.1 

1.3 

5.6 

10.8 

6.6 

Material :  Ove 
B  (cm)  a  (era) 

2.5  5.6 
1.9  5.6 
1.3  S . 6 


219;  AK* 26 . 4 

MPa/m 

Kj  (Mpa/iii) 

K2  (MPa/iii) 

8.2 

3.1 

15.6 

8.8 

11.8 

5.1 

No. 


TABLE  5 
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THE  INFLUENCE  OF  FATIGUE  CRACK 
SURFACE  ROUGHNESS  ON  ACOUSTIC  WAVE  TRANSMISSION 
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The  University  of  Texas  at  Austin 
Austin,  Texas  78712 


0.  Buck 
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Thousand  Oaks,  California 


Abstract 

This  paper  is  concerned  with  t"he  measurements  that  simulate 
the  amplitude  of  acoustic  waves  transmitted  through  a  simulated 
crack  interface  as  the  congressional  load  on  the  crack  is  varied 
and  as  the  roughness  of  the  crack  is  varied.  The  samples 
either  have  specific  roughness  mechanically  produced  or  are  cut 
from  compact  tension  specimen  run  at. different  stress  intensity 
ranges  and  in  different  environments  that  represent  actual 
fatigue  crack  surfaces.  Using  an  electro-hydraulic  testing 
machine,  two  aluminum  blocks  are  loaded  in  compression 
until  the  maximum  acoustic  transmission  is  ooserved.  It  is 
observed  that  the  amplitude  of  acoustic  wave  transmitted  through 
the  Interface  reaches  the  maximum  at  different  stress  levels 
depending  on  the  inter  fa c ia 1  surface  roughness.  Additional 
variables  of  the  study  are  the  frequency  of  the  acoustic  waves, 
and  the  relative  hum;Jity  both  during  the  fatigue  test  and  during 
the  acoustic  evaluation.  A  model  is  described  showing  the  cor¬ 
relations  between  load,  surface  roughness,  acoustic  frequency, 
and  environment  on  the  transmitted  energy.  The  results  are  then 
related  back  to  the  extent  of  crack  closure  observed  acoustically 
on  a  fatigue  cracked  specimen. 


A.  INTRODUCTION 


This  paper  is  concerned  with  the  measure¬ 
ment  of  the  attenuation  of  longitudinal 
acoustic  waves  as  they  are  transmit¬ 
ted  through  a  simulated  fatigue  crack 
interface  in  A1  2219  and  A1  2024  alloys. 
Additionally,  combining  the  acoustic 
data  with  an  elastic  finite  element 


analysis  that  defines  the  stress  dis¬ 
tribution  across  the  interface  it  is 
shown  that  it  is,  in  principle,  pos¬ 
sible  to  determine  the  actual  extent 
of  crack  closure  from  acoustic  data 
taken  during  the  fatigue  cycle.  The 
variables  of  this  study  are  the  ap- 
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plied  load,  surface  roughness  of  the 
simulated  fatigue  crack  interface,  and 
frequency  of  the  acoustic  waves. 


The  machining  led  to  a  much  coarser 
roughness  than  the  grit  blasting. 


B.  SPECIMEN  PREPARATION 

The  aluminum  blocks  were  cut  from  plate 
stock  and  parallel  ground  to  minimize 
any  effects  of  rocking  action  between 
the  surfaces.  After  being  ground,  the 
blocks  were  metal lurgical ly  polished 
on  the  surfaces  that  were  to  serve  as 
the  simulated  interface.  Next,  the  pol¬ 
ished  surfaces  were  grit  blasted  with 
varying  grits  to  produce  the  desired 
variations  in  surface  roughness.  The 
400  microinch(RMS)  sample  was  cut  from 
an  actual  fatigue  crack  grown  in  a 
compact  tension  specimen.  The  authors 
feel  grit  blasted  surfaces  closely  ap¬ 
proximate  the  nature  of  a  real  fatigue 
crack  surface  where  tearing  and  hammer¬ 
ing  have  occurred,  thus  preventing  any 
matching  of  the  surfaces  on  a  micro¬ 
scopic  scale  (Fig.l). 


DEFINED  RMS 
SURFACES 


A*  =  9.6  cm* 


a 


DISCRETE  CONTACT 
POINTS 


f 

a 


FIGURE  1.  THE  SIMULATED  FATIGUE  CRACK 


The  A1  2024  specimens  had  the  model  sur¬ 
faces  machined  to  achieve  the  desired 
roughnesses,  with  the  exception  of  the 
4  microinch  sample  which  was  grit  blasted. 


C.  EXPERIMENTAL  TECHNIQUE 

The  blocks  were  placed  in  an  MTS  closed- 
loop  electrohydraul ic  testing  machine  and 
the  upper  and  lower  mandrels  were  slot¬ 
ted  so  the  longitudinal  wave  acoustic 
transducers  could  be  placed  directly  on 
the  blocks.  The  interface  between  the 
transducers  and  the  blocks  was  of  sil¬ 
icone  based  qrease.  Compression  springs 
were  used  t'  hold  the  transducers  firmly 
against  the  i.ocks.  The  acoustic  signal 
was  generated  by  a  standard  pulse  oscil¬ 
lator,  travelled  from  the  input  trans¬ 
ducer  through  the  simulated  fatigue  crack 
to  the  output  transducer,  then  to  an  AC-DC 
converter,  and  finally  to  an  X-Y  plotter 
where  it  appeared  on  the  ordinate  axis. 

The  load  signal  was  taken  directly  from 
the  MTS  controller  and  appeared  on  the 
abscissa  axis.  Thus  the  outputs  of  the 
experiment  could  be  directly  plotted 
while  the  system  was  cycling  (at  0.03  Hz.) 
Typically  about  fifteen  cycles  were  re¬ 
quired  for  the  output  signal  to  stabilize 
and  the  data  reported  is  from  after  stab¬ 
ilization  of  the  signal  had  occurred. 

The  loading  cycles  were  from  zero  com- 
pressional  load  to  a  compressive  load  of 
approximately  30  percent  of  the  yield 
strength  of  the  material,  based  on  the 
apparent  area  of  contact.  To  seperate 
out  the  transducer  coupling  effect  a 
block  with  no  interface  was  used. 

0.  EXPERIMENTAL  RESULTS 

In  figure  2  the  appearance  of  a  typical 
fatigue  crack  is  shown.  Note  that  part 
of  the  crack  is  open  and  part  closed 
consistent  with  the  fatigue  crack  closure 
model  (1.2)  Typically  acoustic  tech¬ 
niques  are  used  to  discover  cracks  such 
as  this.  Previously  it  has  been  impos¬ 
sible  to  determine  the  size  of  a  fatigue 
crack  from  acoustic  data  because  no  in¬ 
formation  was  available  to  determine  the 
amount  of  acoustic  wave  transmitted'31 
through  the  closed  portion  of  the  crack 
with  crack  topography  as  a  variable. 
Figures  3  and  4  show  the  dramatic  effect 
of  interfacial  surface  roughness  on  the 
amplitude  of  acoustic  wave  transmitted 
through  the  interface  as  a  function  of 
load.  Significantly,  the  difference  in 
transmission  between  the  grit  blasted 
and  machined  surfaces  for  equivalent 
loading  should  also  be  noted.  An  in¬ 
teresting  aspect  of  the  results  from 
the  grit  blasted  specimens  is  what  appears 
to  be  a  saturation  effect  (Fig. 5),  how¬ 
ever,  the  authors  feel  further  invest¬ 
igation  will  be  necessary  to  verify  this 
finding.  Another  variable  of  this  study 
was  frequency  of  the  acoustic  waves. 
Comparison  of  data  taken  for  two  dif- 


ferent  model  cracks  shows  essentially 
no  difference(Fig.  6). 
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result  reported  previously  by  Lohanlc 
and  Llpschul tz  '  '  • 

E.  ELASTIC  FINITE  ELEMENT  ANALYSIS 


FIGURE  2.  A  TYPICAL  FATIGE  CRACK 

Addlonally,  note  the  similarity  in  the 
4  microinch  curves  in  figures  3  and  4. 


The  authors  will  now  give  a  brief  de¬ 
scription  of  the  computer  model  for 
crack  closure  using  a  finite.element 
program  developed  by  Sewell'  .  Assume 
that  a  fatigue  crack  growth  experiment 
is  carried  out  on  a  compact  tension 
specimen.  As  the  crack  grows,  it  will 
leave  behind  an  approximately  uniform 
wake  of  residual  strain,  that  is,  the 
material  in  the  wake  of  the  crack  tip 
will  effectively  become  thicker  because 
of  the  surface  roughness  and  volumn 
change  associated  with  the  high  tensile 
strains  near  the  crack  tip  before  seper 
ation.  Jbis  increase  in  thickness  is 
modeled  6)as  a  uniform  coating  of  new 
material  on  the  crack  surfaces.  The 
normal  stress  distribution  across  the 
crack  interface  is  then  established  as 
a  function  of  distance  from  the  crack 
tip.  The  stress  distribution  obtained 
by  the  analysis  is  then  available  to 
interpret  the  transmitted  acoustic  data 
from  a  fatigue  crackl2). 


FIGURE  3.  ACOUSTIC  WAVE  AMPLITUDE  VS.  APPLIED  STRESS  (GRIT  BLASTED  SPECIMENS) 


The  acoustic  wave  frequencies  used  here 
were  2.25  MHz  for  the  former  and  1.0  MHz 
for  the  latter.  Thus  the  authors  feel 
that  there  is  little  frequency  effect 
in  the  1.0  MHz  to  5.0  MHz  regime,  a 


F.  DISCUSSION  AND  CONCLUSION 


The  authors  have  demonstrated  that  in 
order  to  evaluate  the  true  extent  of 
crack  closure  distance  as  a  function 


FIGURE  A.  ACOUSTIC  WAVE  AMPLITUDE  VS.  APPLIED  STRESS  (MACHINED  SPECIMENS) 
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FIGURE  5.  THE  SATURATION  EFFECT 


of  load  it  is  necessary  to  consider  the 
compressive  stress  on  the  fracture  sur¬ 
face,  an d  the  transmission  factor 
which  is  a  function  of  the  stress  and 
surface  roughness.  Usis>g,the  elastic 
finite  element  analysis'  ,  the  extent 
of  closure  and  the  closure  stress  have 
been  determined  allowing  the  calcula¬ 
tion  of  the  acoustic  wave  transmission 
across  the  interface  which  s  given  by; 


where  T^j  =  total  acoustic  transmission 
across  the  interface 
a=  true  crack  length 


:rir< 


x ) RMS ) d  x 


•vV*? 


x»  distance  along  crack  from 
load  line  to  a. 

and, 

T(e(x)RMS)’  the  transmission  factor 
as  a  function  of  stress 
and  surface  roughness 

Therefore  it  is.  In  principle,  possible 
to  predict  the  amount  of  transmission 
due  to  crack  closure  from  the  roughness- 
stress  data  presented  here  combined  with 
the  use  of  the  elastic  finite  element 
model.  It  should  be  noted  that  the  above 
analysis  is  for  the  instantaneous  measure¬ 
ment  of  crack  lengths.  For  NOT  purposes 
the  non-elastic  behavior  of  the  fatigue 
crack‘7l  must  be  taken  into  consideration. 


FIGURE  6.  THE  FREQUENCY  EFFECT 


Another  application  of  the  results  of 
this  study  is  to  the  choice  of  load 
required  to  acoustically  couple  a  re¬ 
mote  sensing  transducer  through  a  sens¬ 
ing  rod  to  the  structure  of  Interest. 

In  this  case  the  RMS  should  be  minimized 
to  reduce  the  required  loading. 
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MODEL  FOR  ULTRASONIC  MEASUREMENT  OF  TRUE 
CRACK  LENGTH  WHEN  CRACK  CLOSURE  EXISTS 
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MODEL 

The  amplitude  of  acoustic  wave  transmission  through  the  closed 
portion  of  a  fatigue  crack  is  a  function  of  two  variables.  These  vari¬ 
ables  are  the  roughness  of  the  fatigue  crack  interface  and  the  compres¬ 
sive  stress  distribution  from  the  point  of  crack  closure  to  the  crack  tip. 
The  model  to  be  developed  will  encompass  both  of  these  variables. 

In  general,  the  curves  described  by  the  plots  of  acoustic  wave 
transmission  versus  applied  stress  (Fig.  4),  are  roughly  approximated  by 
the  form: 

A  •  'Wd-e*'”)  -1- 

where  A  =  the  amplitude  of  the  transmitted  acoustic  wave  normalized  to  the 
maximum  transmitted  wave 

Amax  =  maximum  transmitted  wave 

n  =  a  variable  to  be  determined 

and  o  *  the  stress  across  the  interface  which  is  always  compressive. 

For  the  roughnesses  produced  in  the  compact  tension  specimens, 


it  was  found  that  the  curve  described  by: 
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A  =  Amax(l-e+'028a)  -2- 

closely  matches  the  acoustic  wave  amplitude  versus  applied  stress  curve 
previously  measured  (Fig.  10).  Here  the  units  of  A  are  arbitrary  with  a 
maximum  of  very  nearly  and  a,  the  applied  stress,  is  in  units  of 
MN/m2. 

If  it  is  assumed  that  the  stress  distribution  is  linear  and 
increases  from  zero  at  the  point  of  closure  to  the  yield  strength  of  the 
material  at  the  crack  tip  then  it  is  of  the  form: 

a  =  mx  -3- 

where  m  is  a  constant  dependent  on  x,  the  crack  closure  distance,  and  x 
varies  from  zero  at  the  crack  closure  point  to  xT,  the  distance  from  the 
closure  point  to  the  crack  tip. 

The  expression: 

.  .  . ,  +mnx , 

A  ~  Amax( 1-e  )  -4- 

describes  the  predicted  transducer  output  curve.  Next,  the  values  of  m 
for  the  two  specimens  must  be  determined.  For  specimen  one: 
at  x  =  0,  o  =  0 

at  x  =  xT  =  .018m,  o  =  -380MN/m2 

°°°  m  =  *  =  Tuff  =  ”21 , 100MN/m3 

For  specimen  two: 

at  x  =  0,  o=0 

at  x  =  xT  =  .005m,  a  =  -380MN/m2 
o°o  m2  =  f  =  ^  =  -76,OOOMN/m3 

where  the  location  of  x  =  0  and  x  *  xT  have  been  determined  experimentally. 
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Thus  the  two  equations  describing  the  normalized  predicted  acoustic  wave 
output  for  the  two  specimens  are: 


A,  * 

-5- 

and 

A,  *  AInax(l-e4OTnx)  -  -?2'50*) 

-6- 

In  order  to  determine  the  total  normalized  acoustic  wave 
across  the  interface,  Tai»  the  curves  must  be  integrated 
x  =  xT,  or 


x=x 

Tai  =  Amax  f~ 
x=0 


T 

1-e 


+mnx 


dx  =  Amax 


+mnx 

e 


transmission 
from  x  =  0  to 


Evaluation  of  the  integral  for  the  two  curves  yields: 


TAh“ 

TAI2  =  -OOABAmax 

and  the  units  for  Taii  and  TA12  will  be  distance  times  the  maximum  wave 
transmission. 

In  order  to  determine  the  perc  :ntage  of  acoustic  wave  transmis¬ 
sion  across  the  interface,  TAI1  and  TAI2  must  be  divided  by  the  integral 
of  the  acoustic  wave  if  all  transmission  began  at  the  closure  point,  or: 


x=xT 

^AT  =  f  Amax^X 
x=0 


-8- 


where  T  is  the  total  acoustic  wave  transmission  assuming  all  transmis¬ 
sion  began  at  the  closure  point  and  A  represents  the  maximum  trans¬ 
mission  across  the  interface.  Thus 

x=xT  x=xT 

^AIl  “  f  ^max  dx"Atnax  !  “  •  018Atnax  -9- 

x=0  x=0 


AI2 


x=xT  x=xT 

/  Amax  dx=Amaxx  I  =  -005\ax 


x=0 


x=0 
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where  the  units  will  again  be  distance  times  the  maximum  acoustic  wave 
transmission.  Now  the  transmission  percentage,  T  ,  can  be  obtained  and 
is  equal  to 


'AP 


Therefore 


Iap2  = 


. 0l63Amax 

•01 8U 
.0045^ 

■OOBAjnax 


=  .90  =  90% 


=  .90  =  90% 


and  the  units  of  the  numerator  and  denominator  cancel  each  other.  Thus 
10%  of  the  acoustic  wave  transmission  is  attenuated  by  the  low  stress 
portion  of  the  crack. 

In  the  preceeding  section  of  the  model  development  it  was 
assumed  that  the  stress  distribution  from  the  crack  closure  point  to  the 
crack  tip  was  linear.  However,  it  is  better  to  model  the  stress  distri¬ 
bution  as 

o  =  mx2  -11- 

where  m  is  a  constant  and  dependent  upon  the  degress  of  crack  closure. 

This  stress  distribution  approximates  the  results  of  an  elastic  finite 
element  analysis  in  Appendix  E  or  that  expected  from  beam  theory. 

The  values  of  m  for  the  nonlinear  stress  distribution  model  are  obtained 
in  the  same  manner  as  for  the  linear  model. 


For  specimen  one: 

at  x  =  0,  o=0 


Therefore,  the  expression  for  the  total  normalized  acoustic  wave 
transmission  across  the  interface  for  the  two  specimens  are 


x=.018 

Taii  *  <wj  (l-e‘3-33xl0‘'x2)dx  -16- 

x=0 
x= . 005 

and  Tai2  =  wj  (l-e‘4*31xl05x)dx  -17- 

x=0 

which,  upon  evaluation,  yield 

Taii  =  1*31  x  10  A,nax 

and  _~ 

Tai2  =  3.65  x  10  Ajnax 

The  units  of  TAI1  and  Tai2  will  again  be  distance  times  the  maximum 
acoustic  wave  transmitted.  To  obtain  the  percentage  of  wave  transmission 
across  the  interface,  the  amount  of  transmission  that  would  occur  if  all 
transmission  began  at  the  closure  point  must  be  determined.  This  quantity 
TAt,  is  equal  to  that  calculated  in  equations  8  -  10.  Now  the  trans¬ 
mission  percentage,  TAP,  can  be  determined. 


TAI1 

T  - 

_  1.31xl0"2A 

.73  =  73% 

-18- 

API  T 

‘ah 

.018A 

and 

t  _  TAI1 

^  TAI2 

_  3.65xl0"3Amax  = 

. 005Amax 

.73  =  73% 

-19- 

Therefore,  27%  of  the  acoustic  wave  transmission  is  attenuated 
by  the  low  stress  portion  of  each  crack.  And  it  has  been  shown  that  the 
percentage  of  acoustic  wave  transmission  is  independent  of  the  crack 
closure  distance.  Note  that  the  percentage  transmission  obtained  is  only 
for  the  conditions  of  o  =  mnx2  stress  distribution  and  a  transmission  co¬ 
efficient,  A,  independent  of  surface  roughness  over  the  range  of  interest. 


It  was  demonstrated  in  Appendix  c  that 
acoustic  wave  transmission  through  a  fatigue  crack  interface  is  virtually 
independent  of  roughness  in  the  4.8  pmeter  to  10  pmeter  rms  range.  In  the 
development  of  the  model  in  the  preceeding  section,  it  was  shown  that  the 
percentage  of  acoustic  wave  transmission  through  a  fatigue  crack  is  inde¬ 
pendent  of  closure  distance  for  a  fixed  stress  profile.  This  independence 
of  closure  distance  applies  to  both  the  linear  and  quadratic  models  of  the 
stress  distribution  normal  to  the  surface  along  the  closure  distance.  In 
the  prediction  of  true  fatigue  crack  size  the  model  can  be  used  in  several 
ways. 

First  of  all,  it  must  be  noted  that  when  the  crack  analysis  uses 
measurement  of  reflected  waves  rather  than  transmitted  waves,  the  pre¬ 
dicted  transmission  quantity  must  be  subtracted  from  one  to  give  the 
reflected  wave  quantity. 

The  two  most  general  cases  of  crack  interrogation  occur  when  the 
transducers  are  large  relative  to  the  closure  distance  and  when  point 
transducers  are  used  (See  Fig.  2). 

In  the  use  of  transducers  large  relative  to  the  closure  dis¬ 
tance  two  different  techniques  may  be  used  to  account  for  the  crack  clos¬ 
ure  and  to  predict  the  crack  size.  The  first  of  the  two  techniques 
requires  the  use  of  the  predicted  percentage  of  acoustic  wave  transmis¬ 
sion.  Since  the  transducer  is  large  relative  to  the  closure  distance  the 
output  can  be  considered  the  integral  of  the  acoustic  wave  transmitted 


x=o 


LARGE  TRANSDUCER  RELATIVE  TO  CLOSURE 


— 

0 
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FIGURE  2.  THREE  TRANSDUCER  SIZES  RELATIVE  TO  CRACK 
CLOSURE  DISTANCE 


When  point  transducers  are  used  to  examine  a  fatigue  crack,  the 
procedure  is  essentially  the  same  as  that  described  for  the  large  trans¬ 
ducers.  Instead  of  locating  the  point  where  the  change  in  output  becomes 
linear,  the  point  where  transmission  first  reaches  constant  output  corre¬ 
sponds  to  the  point,  xc,  where  the  stresses  are  high  enough  to  make  the 
crack  appear  transparent  to  the  acoustic  waves.  Now  the  distance  from 
the  closure  point  to  the  crack  tip  may  be  determined  just  as  in  the  pro¬ 
ceeding  section  for  large  transducers. 

When  transducers  are  used  in  fatigue  crack  analysis  that  are 
smaller  than  the  closure  distance,  but  larger  :han  point  transducers  and 
have  a  non-uniform  front,  the  evaluation  becomes  more  difficult.  Now,  not 
only  the  acoustic  wave  transmission  as  a  function  of  stress  can  be  deter¬ 
mined,  but  the  geometry  of  the  transducer  itself  may  need  to  be  considered. 
There  are  two  techniques  that  may  be  used  with  transducers  in  this  size 
range.  The  first  technique  is  similar  to  that  developed  for  large  trans¬ 
ducers  when  the  predicted  percentage  of  transmission  is  used.  Since  the 
transducers  are  smaller  than  the  crack  closure  distance,  however,  as  the 
transducers  are  moved  along  the  crack  the  outputs  must  be  summed.  This 
summation  should  provide  the  equivalent  of  the  integral  of  the  acoustic 
wave  output  that  would  be  observed  if  large  transducers  were  used  when 
the  geometric  effect  is  taken  into  account.  The  technique  now  is  that 
developed  earlier  for  the  large  transducers. 

The  second  technique  for  use  with  these  transducers  is  essenti¬ 
ally  the  same  as  the  method  described  for  use  with  the  point  transducers. 
The  point  where  transmission  begins  still  corresponds  to  the  crack  closure 
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through  the  crack.  The  expected  output  of  the  transducer  is  shown  in 
Fig.  3,  At  some  point  along  the  output  curve  it  becomes  linear.  This 
linear  portion  of  the  curve  corresponds  to  the  increase  in  transducer 
area  parallel  to  the  crack,  past  the  point  where  the  crack  appears  trans¬ 
parent  to  the  acoustic  waves.  If  a  line  is  constructed  beginning  at 
x  *  0,  the  closure  point,  with  the  same  slope  as  the  linear  portion  of 
the  output  curve,  this  new  curve  represents  the  integral  of  acoustic  wave 
transmission  if  complete  transmission  occurred  at  the  closure  point  (i.e., 
it  is  the  integral  of  A^x).  Thus,  from  the  model,  when  the  output 
.curve  is  the  predicted  percentage  of  the  integral  of  the  A^x  curve  the 
crack  tip  has  been  located.  An  example  of  this  is  shown  in  Fig.  12. 

The  second  technique  for  use  with  the  large  transducers  requires 
the  location  of  the  point  along  the  output  curve  where  the  curve  first 
becomes  linear.  This  point  corresponds  to  the  point  along  the  crack 
where  transmission  is  complete.  Experimentally,  by  knowing  the  point 
where  transmission  begins  (x  =  0),  the  point  (xc)  and  the  stress  (,24ay) 
where  transmission  is  95%  complete  and  assuming  the  stress  at  the  crack 
tip  (xT)  is  equal  to  the  yield  strength  of  the  material,  the  crack  length 
can  be  determined.  To  accomplish  this,  the  assumed  stress  distribution 
o  =  nx2 

is  used  where  x  is  the  distance  from  the  closure  point  to  the  crack  tip. 
The  value  of  n  is  then  determined  by  the  point  where  transmission  is  95% 
complete.  Then  the  crack  length  can  be  determined  from  the  stress  of  <?y 
at  the  crack  tip. 


•p. 


X=0  Is  the  crack  closure  point 

X=XC  is  the  point  where  acoustic  wave  transmission  is  complete 
X=X.  is  the  location  of  the  predicted  crack  tip 
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point  where  the  compressive  stress  is  zero.  To  determine  the  point  where 
acoustic  transmission  is  complete  the  location  of  the  trailing  edge  of 
the  transducer  must  be  known.  When  the  acoustic  transmission  reaches  a 
constant  value,  the  trailing  edge  of  the  transducer  is  at  the  point,  xc, 
where  transmission  is  complete.  Therefore,  as  in  the  case  of  the  point 
transducers,  enough  information  is  known  about  the  crack  to  predict  the 
true  crack  length. 

A  model  has  now  been  developed  and  described  which  allows  the 
prediction  of  fatigue  crack  size  in  aluminum  alloys  for  three  different 
size  transducers  relative  to  the  crack  closure  distance.  The  model  should 
also  be  applicable  to  fatigue  crack  analysis  in  other  materials  when  the 
required  data  are  obtained. 


CONCLUSION 


In  this  research,  a  model  has  been  developed  and  described  which 
allows  the  prediction  of  true  fatigue  crack  size  based  on  acoustic  trans¬ 
mission  analysis  in  aluminum  alloys.  Additionally,  it  has  been  shown  that 
the  acoustic  wave  transmission  through  a  fatigue  crack  interface  as  a 
function  of  applied  stress  is  very  dependent  on  the  interfacial  surface 
roughness  when  the  roughness  is  below  approximately  4.8  umeters  rms,  and 
independent  of  roughness  in  the  4.8  umeter  to  10  umeter  rms  range.  The 
10  umeter  rms  roughness  was  the  maximum  roughness  tested  and  corresponds 
to  measured  values  from  samples  fatigued  with  a  constant  stress  intensity 
range. 

This  model  should  be  applicable  to  fatigue  crack  analysis  in 
other  materials  such  as  steel  and  titanium  when  the  required  data  are 
obtained.  The  model  predicts  the  crack  length  independent  of  relaxation 
of  the  crack  as  long  as  the  general  stress  dependence  across  the  crack 
interface  remains  fixed  and  the  crack  surfaces  are  not  totally  in  contact. 
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APPENDIX  g-1 


Evaluation  of  the  integral 


A  =  /(i-e'ax2)dx 
x=0 


The  curve  of  1-e"  describes  a  normal  distribution  and  must  be  solved 


by  approximation. 


x  -ax2 

A  =  ( (1-e  x  )dx 
x=0 


Putting  in  the  finite  value  of  x  for  the  upp;r  limit  of  integration,  let 


the  upper  limit  be  n 


A  =  n-  j  (l-e~ax2)dx 
x=0 


Let  y  =  /5x  then  y2  =  ax2 


and  dy  =  /a  dx 


A  —  / 


Let  z  =  /2y  then  z  =  2y2 

and  dz  =  JT  dy 

n  /T/T 


■  w  /■ 


e  ‘  '2  dz 

JT 
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nja-  JT 

f  e'z2/2dz 

0 


A  =  n 


.  i_  JE 

J2a  1 


0 


A  =  n 


if  JTa 


>  V  then  a 
D 


50 


(1)  if  a  >  50  then  A  -  n  -  1-& 

2  /a 

(2)  ifa<5Q  then  A  =  n  -  N(n  STa) 

and  N  must  be  found  in  a  *.ible  for  normal  distributions. 

In  the  development  of  the  model,  a  is  always  greater  than  n1 2/50  so 
the  approximation  (1)  is  used. 
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Communications 

Oxygen  Transport  During  Fatigue 
Crack  Growth 

J.  W.  SWANSON  AND  H,  L.  MARCUS 

It  is  generally  observed  that  fatigue  crack  growth 
rate  in  an  active  gaseous  environment  is  higher  than 
growth  rate  in  a  baseline  vacuum  environment.1’2  It  is 
not  clear  at  this  point  whether  this  is  a  result  of  inter¬ 
action  at  the  crack  tip  in  the  form  of  a  cyclic  stress 
corrosion  or  if  the  gaseous  species  are  transported 
into  the  material  and  internal  interaction  with  the  de¬ 
formation  behavior  occurs.  In  the  present  work,  com¬ 
parison  of  fatigue  crack  propagation  rates  in  vacuum 
and  oxygen  environments  was  made  for  Monel  404,* 

•Rcgistetcd  trademark.  The  international  Nickel  Company,  Inc. 

commercially  pure  titanium,  and  aluminum  2219-T87 
and  7075-T651  alloys  (see  Table  I  for  compositions). 
Following  the  fatigue  experiments,  samples  of  the  fa¬ 
tigue  surfaces  were  analyzed  with  Auger  Electron 
Spectroscopy  (AES)  and  Secondary  Ion  Mass  Spectros¬ 
copy  (SIMS)  to  determine  if  oxygen  had  been  trans¬ 
ported  into  the  material  during  crack  propagation.  The 
AES  measurements  were  used  due  to  the  better  spatial 
resolution  and  consistent  surface  chemical  sensitivity. 
The  SIMS  allowed  isotope  discrimination  during  the 
sputter  profiling  of  oxygen  into  the  bulk. 

Compact  tension  specimens  were  fatigue  cycled 
through  metal  bellows  in  an  ion-pumped  vacuum  cham¬ 
ber.  Al  a  pressure  of  1.3  pPa  (10"“  Torr),  the  crack 
was  grown  0.8  cm  with  load  shedding  accomplished 
each  0.16  cm  of  growth  to  maintain  me  AE  relatively 
constant  at  the  value  given  in  Table  I.  Growth  rate 
was  thus  established  for  a  vacuum  environment.  Under 
constant  AE  conditions,  crack  growth  rate  should  be 
constant  except  fov  environmental  effects,  if  present 
The  isotope  oxygen-18  of  99.9  pet  purity  was  intro¬ 
duced  into  the  chamber,  backfilling  to  a  pressure  of  10 
kPa  (1  kPa  for  aluminum  2219  due  to  a  limited  supply 
of  the  isotope).  The  crack  was  extended  another  0.8 
cm  at  the  fixed  A E  and  growth  rate  in  the  oxygen- 18 
determined.  The  specimens  were  cut  and  transferred 
to  the  AES/SIMS  system.  Compositional  depth  pro¬ 
files  of  oxygen  into  the  fatigue  crack  surfaces  were 
then  obtained  by  use  of  inert  ion  sputter-etching,  com¬ 
bined  with  AES  and  SIMS.  A  more  detailed  description 
ol  the  apparatus  and  procedure  can  be  found  in  Ref.  3. 

The  surface  formed  in  vacuum  had  been  exposed  to 
oxygen- 18  longer  thanthe  surface  fatigued  in  the  oxygen- 
18.  since  it  was  formed  prior  to  the  introduction  of  the  oxy¬ 
gen-  18  K  is  expected  that  the  profile  from  the  frac- 
t  if  surface  formed  in  vacuum  represents  simply  oxi- 
i,,n  ami  diffusion,  while  the  fracture  surface  pro- 
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file  from  the  sample  fatigued  in  oxygen- 18  represents 
oxidation,  diffusion,  and  any  additional  oxygen  trans¬ 
port  mechanism  occurring  during  the  fatigue  crack 
propagation.  Oxidation  and  diffusion  that  occurred 
during  transfer  from  the  fatigue  chamber  to  the  analy¬ 
sis  chamber  is  in  predominantly  oxygen-16. 

Results  of  the  fatigue  crack  growth  experiments 
showed  an  increased  crack  growth  rate  in  the  oxygen 
environment  by  a  factor  of  5  and  2  for  Monel  and  ti¬ 
tanium,  respectively,  with  little  or  no  increase  ap¬ 
parent  for  the  aluminum  alloys.  The  penetration  depth 
of  the  oxygen  in  each  case  as  measured  by  AES  and 
SIMS  is  shown  in  Figs.  1  through  8.  It  was  greater  in 


Sputtering  Time  (Min) 

Fig.  1— AES  oxygen  profiles  from  Monel  404  fatigue  crack  sur¬ 
faces  formed  in  vacuum  and  oxygen-18  environments. 
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Fig.  2— SIMS  oxygcn-18  profiles  from  Monel  40 1  fatigue  crack 
surfaces  formed  in  vacuum  and  oxygen-18  environments. 


Table  I.  Nominal  Malarial  Compotitiom  and  Constant  Sutu 
Intensity  Ranga  Usad  in  Fatigua  Crack  Growth  Studiat 


Material 

Composition  (Pet) 

AA(MPaZTh) 

Monel  404 

52-57  Ni.  0  5  max.  Fe, 

0.1  max.  Mn,  0.1 5  max. 

C.  0  024  max.  S,  0.1 
max.  Si.  0.05  max.  Al. 
bal  Cu 

JJ 

UfanJum 

(Yminicrclally  pure, 

ASTM  B  265  (Cr  2). 

99  2  Ti 

7  7 

Aluminum  707  *y  ThM 

5  6  In.  2  5  Mg.  O.M'u. 

0.3  IT.  bal,  Al 

11 

Aluminum  2Jl‘>-T*7 

6.1  Cu.  0.3  Mn.  0.06  Ti. 

0  1  V.  0.2  Zr.  bal  Ai 

11 

" .%  \ 
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Fig.  3—AES  oxygen  profiles  from  titanium  fatigue  crack  sur¬ 
faces  formed  in  vacuum  and  oxygen-18  environments. 
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Fig.  4 — SIMS  oxygen-18  profiles  from  titanium  fatigue  crack 
surfaces  formed  in  vacuum  and  oxygen-18  environments. 


Fig.  5— AES  oxygen  profiles  from  2219  aluminum  fatigue  crack 
surfaces  formed  in  vacuum  and  oxygen- 18  environments. 

all  cases  when  the  crack  was  grown  in  oxygen- 18. 

There  is  consistency  between  AES  and  SIMS  in  indi¬ 
cating  enhanced  penetration.  Limited  sputter  rate 
calibration  prevents  quantification  of  transport  depths, 
however  using  an  estimated  50  to  100A  for  the  natural 
oxide  thickness  gives  an  estimate  of  the  enhanced 
transport  of  100  to  600A. 

Samples  of  metallographieally  polished  Monel  and 
titanium  were  allowed  to  oxidize  and  were  sputter  pro¬ 
filed  using  the  same  Ion  beam  parameters  as  for  the 
fatigue  samples.  The  profiles  thus  obtained  demon¬ 
strated  point  by  point  agreement  with  the  samples  fa¬ 
tigued  in  vacuum.  This  would  suggest  that  changes  in 
surface  roughness  associated  with  fatigue  crack  growth 
did  not  contribute  significant  error  to  the  measurement 
and  could  not  explain  the  differences  observed  in  the 
vacuum  and  oxygen- 18  runs. 
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Fig.  6— SIMS  oxygcn-18  profiles  from  2219  aluminum  fatigue 
crack  surfaces  formed  :i  vacuum  and  oxygen-18  environments. 
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Fig.  7— AES  oxygen  profiles  from  7075  aluminum  fatigue  crack 
surfaces  formed  In  vacuum  and  oxygen-18  environments. 
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Fig.  8— SIMS  oxygen-18  profiles  from  7075  aluminum  fatigue 
crack  surfaces  formed  in  vacuum  and  oxygen-18  environments. 

From  scanning  electron  micrographs,  no  change  in 
fracture  surface  appearance  was  noted  for  the  alumi¬ 
num  alloys  or  titanium.  The  Monel  exhibited  a  targe 
percentage  of  fracture  of  a  crystallographic  nature 
when  fatigued  in  oxygen,  in  agreement  with  previously 
reported  results.' 

A  model  for  sweep- in  of  interstitial  atoms  that  ori¬ 
ginate  in  the  gas  phase  by  mobile  dislocations  during 
fatigue  crack  propagation  has  been  proposed5’0  and 
shows  the  plausibility  of  such  transport  for  the  small 
hydrogen  atom.  Using  the  bulk  diffusivity  of  oxygen, 
this  model  is  Insufficient  to  explain  the  present  re¬ 
sults.  If  an  enhanced  oxygen  diffusivity  near  the  re¬ 
laxed  high  dislocation  density  surface  is  postulated, 
then  the  dislocation  mechanism  is  plausible. 

The  interpretation  of  the  fatigue  crack  growth  re¬ 
sults  by  the  authors  is  that  gaseous  environmental 
species,  oxygen  in  this  case,  are  transported  into  the 
metal  during  passage  of  a  fatigue  crack,  possibly  by 
mobile  dislocations.  The  interstitial  solute  atoms  In 
increased  concentrations  then  interact  with  the  base 
metal,  significantly  modifying  the  local  plastic  be¬ 
havior  and  fracture  strain  in  the  vicinity  of  the  crack 
tip.  This  interaction  for  oxygen  is  apparently  the 
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strongest  for  Monel  404  and  minimal  for  the  alloyed 
aluminum.  The  oxygen  influence  on  fracture  strain 
could  be  the  cause  for  the  increase  in  growth  rate.  A 
second  possibility  would  be  associated  with  the  crack 
closure  concept’  which  is  concerned  with  the  fact 
that  during  part  of  the  loading  cycle  the  crack  closes 
due  to  a  residual  displacement  pattern  in  the  wake  of 
the  crack.  In  the  case  reported  here  the  observed 
change  in  oxygen  penetration  and  postulated  change  in  , 
plastic  behavior  due  to  the  transport  of  the  oxygen 
could  be  expected  to  result  in  a  smaller  residual  dis-  • 
placement  in  the  wake  of  the  crack,  increasing  the  ef~ 
fective  A K  used  in  the  crack  closure  model,  conse¬ 
quently  increasing  the  growth  rate.  A  great  deal  of 
additional  research  will  be  required  to  establish  the 
validity  of  these  hypotheses. 
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1.  Introduction 


The  purpose  of  this  investigation  was  to  study  the  effect  of  hydrogen,  deuter¬ 
ium,  water  vapor  and  oxygen  on  the  fatigue  crack  growth  of  high  strength  alu¬ 
minum  alloys.  Fatigue  crack  growth  tests  in  a  controlled  environment  were 
followed  by  Secondary  Ion  Mass  Spectroscopy  (SIMS)  and  Auger  Electron 
Spectroscopy  (AES). 

Influence  of  humid  air  on  fatigue  crack  growth  in  high-strength  aluminum 
alloys  is  an  important  practical  problem  well  documented  in  literature  [1-5]. 
The  studies  show  an  increase  of  crack  growth  rate  of  aluminum  alloys  in  water 
vapor.  In  many  cases  it  is  attributed  to  hydrogen  embrittlement.  On  the 
other  hand  gaseous  hydrogen  has  no  effect  on  fatigue  crack  growth  [4], 

Several  models  of  hydrogen  embrittlement  for  various  materials  have  been 
considered  [6-10].  The  sweep-in  mechanism  during  fatigue  crack  growth  pro¬ 
poses  combination  of  diffusion  and  dislocation  sweeping  of  hydrogen  into  the 
plastic  zone.  This  results  in  a  hydrogen  penetration  deeper  than  JBt,  the 
characteristic  depth  for  diffusion,  and  a  consequent  increase  in  hydrogen 
concentration  that  in  some  manner  weakens  the  material  [9,10].  The  hydrogen 
effect  is  even  more  pronounced  in  the  presence  of  various  trapping  sites 
such  as  particle  interfaces,  which  are  made  more  effective  by  the  deformation 
associated  with  the  propagating  crack  [9,10,11]. 

2.  Experimental  Procedure 

The  age  hardenable  aluminum  alloys  used  in  this  study  are  7075-T651  (chemical 
composition  weight  percent:  5.6%  Zn;  2.5%  Mg;  1.6%  Cu;  0.3%  Cr;  and  cy=505  MPa; 
Kjr=25  MPavfii)  and  2219-T851  (chem  comp.:  6.3%  Cu;  0.3%  Mn;  0.18%  Zn;  0.1%  V; 
0.06%  Ti;  and  Oy=350  MPa;  Ktq=35  MPa/m).  The  compact  tension  fracture  mechan¬ 
ics  type  specimens  were  used  in  the  fatigue  crack  growth  experiments  [12], 

The  fatigue  crack  was  first  grown  in  a  1  yPa  vacuum  at  a  constant  AK.  (ranging 
from  10  to  21  MPavid)  and  frequency  (6  or  10  Hz).  The  chamber  was  than  filled 
with  one  of  the  gases  hydrogen,  water  vapor,  oxygen-18,  deuteriurated  water 
vapor  or  deuterium  gas.  Fatigue  crack  growth  was  then  continued  in  the  en¬ 
vironment,  up  to  a  certain  crack  length  and  then  the  specimen  was  fast  frac¬ 
tured.  Crack  propagation  rate  was  monitored  optically.  Specimens  cut  from 
the  fracture  surfaces  were  then  transferred  to  AES/SIMS  vacuum  system.  Hydro¬ 
gen,  deuterium  and  oxygen  18  concentration  profiles  normal  to  the  fracture 
surface  were  investigated  by  inert  ion  sputtering  combined  with  SIMS.  Nega¬ 
tive  SIMS  was  utilized  in  case  of  deuterium  due  to  the  higher  yield  of  D"  over 
D+  ions.  For  oxygen  AES  was  used  in  parallel.  In  addition  deuterium  ion 
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implant  ed  standards  were  prepared  using  a  high  energy  ion  implantation.  These 
were  used  for  calibration  purposes  of  the  SIMS  system  and  as  a  source  of 
measuring  induced  trapping  effects. 

3.  Results  and  Discussion 


Experimental  difficulties  were  encountered  in  measuring  the  atomic  hydrogen 
profile  with  SIMS  in  the  presence  of  the  3xl0'5  torr  argon  sputtering  gas. 
Therefore,  oxygen  18  and  deuterium  as  isotope  markers  in  SIMS  profiling  were 
used. 


Fatigue  crack  propagation  rate  for  dry  oxygen  18  did  not  significantly 
increase  in  comparison  to  vacuum  (1.1  times).  The  ion  sputtering  profiles 
determined  by  AES  and  SIMS  for  2219  and  7075  aluminum  alloys  have  been  pub¬ 
lished  elsewhere  [12,13].  A  significant  enhancement  of  the  thickness  of  the 
oxygen  for  the  region  fatigued  in  oxygen  18  by  a  factor  of  about  2  to  that 
fatigued  in  vacuum  and  subsequently  exposed  to  oxygen  18  occurred  in  both 
cases.  Dry  hydrogen  or  deuterium  also  show  no  change  in  fatigue  crack  growth 
rate  when  compared  to  vacuum. 

Water  vapor  (H?0  or  D2O)  strongly  influences  the  fatigue  crack  growth  rate. 
The  ratio  of  crack  growth  rate  in  water  vapor  to  that  in  vacuum  varied  from 
1.8  to  3.6  depending  on  the  AK  value.  The  effect  was  more  pronounced  for 
intermediate  AK  values. 


The  SIMS  profile  of  deuterium  in  the  fracture  surface  for  the  specimen 
tested  in  vacuum  represents  deuterium  for  which  diffusion  occurred  when  this 
surface  was  exposed  to  D2O  atmosphere  subsequent  to  the  fatigue  crack  growth 
(Fig. 1 ,  curve  1).  The  surface  profile  for  the  segment  cracked  in  the  D2Q 
vapor  represents  diffusion  and  any  additional  deuterium  transport  mechanism 
occurring  during  the  fatigue  cycling  process  (Fig.l,  curve  2).  Deeper  pene¬ 
tration  of  deuterium  for  the  portion  of  specimen  where  fatigue  crack  growth 
was  in  D2O  vapor  to  that  in  vacuum  was  observed  (Fig.l).  This  is  in  quali¬ 
tative  agreement  with  dislocation  sweeping  mechanism  [9,10]. 


Fio.l  SIMS  profiles  of  D"  in  Al-  Fig. 2  D"  SIMS  profiles  of  deuterium 
2219  T851  alloy  on  fractured  ion  implanted  Mg-base  Mg-Zn  alloy, 

surfaces  and  dummy  specimen. 


The  deuterium  background  level  in  the  AES/SIMS  system  was  obtained  using 
the  dummy  specimen.  This  specimen  was  exposed  to  D2O  environment  for  the  same 
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period  of  time  as  the  others  but  was  not  fatigued  and'  only  a  prior  existing 
oxidized  surface  was  exposed.  The  deuterium  SIMS  profile  from  the  dummy  spec¬ 
imen  is  initially  only  about  10%  of  the  signal  obtained  from  the  D?0  and 
vacuum  fractured  specimens  and  disappears  within  3  tc  4  minutes  of“sputtering 
time  (Fig.l,  curve  3).  This  shows  that  deuterium  contamination  from  the 
AES/SIMS  vacuum  is  very  low. 

Preliminary  attempts  at  calibrating  the  SIMS  system  were  done  on  the  Mg- 
base  Mg-Zn  alloy  ion  implanted  specimens.  •  To  simulate  the  fatigue  fracture 
surface  conditions,  the  specimens  were  ion  bombarded  with  neon  at  energies 
130  keV  and  150  keV.  By  such  multi-implantation,  radiation  damage  was  deep 
and  uniform.  Then  a  known  dose  of  deuterium  was  implanted  into  the  specimen 
within  the  damaged  region.-  The  profiles  of  deuterium  determined  by  SIMS  are 
shown  in  Fig. 2.  Fig. 2  provides  two  other  profiles  taken  two  and  four  days 
later.  These  experiments  will  be  repeated  in  the  near  future  using  aluminum 
alloys  of  our  interest.  Hydrogen  implanted  into  the  single  crystal  of  alu¬ 
minum  has  been  reported  to  diffuse  out  within  half  an  hour  [15].  Therefore, 
presence  of  trapping  sites  is  required  for  retarding  the  diffusion  of  hydrogen 
and/or  deuterium  atoms  [11,14].  During  the  deformation  associated  with  the 
plastic  zone  of  propagating  crack  in  fatigue  fractured  specimens  large  amounts 
of  various  types  of  trapping  sites  are  created.  Using  radiation  damage  we 
hope  to  simulate  the  destruction  of  the  material  of  the  character  created  in 
fatigue  fracture.  This  could  then  be  used  to  evaluate  the  hydrogen  present 
and  its  transport  mechanism.  . 

4.  Conclusions 

(1)  SIMS  can  be  applied  to  study  the  influence  of  water  vapor  on  the  fatigue 
crack  growth  of  aluminum  alloys.  (2)  Dry  hydrogen  (or  deuterium)  does  not 
affect  fatigue  crack  growth  in  2219  and  7075  aluminum  alloys.  (3)  Oxygen 
penetrates  deeper  in  surfaces  formed  during  fatigue  than  vacuum  fatigued  and 
oxygen  exposed  surfaces.  It  does  not  play  a  significant  role  in  fatigue 
crack  growth.  (4)  Hydrogen  (or  deuterium)  from  water  vapor  influences  the 
fatigue  crack  growth  rate  and  increases  penetration  of  hydrogen  and/or  deuter¬ 
ium  into  the  material.  (5)  Ion  implanted  specimens  can  be  used  for  the  SIMS 
calibration  and  for  the  trapping  analysis. 
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Introduction 

The  aim  of  this  paper  is  to  report  on  a  deuterium  in 
aluminum  trapping  inversion  experiment.  Deuterium  ions, 
deposited  at  a  sharply  defined  depth  by  cyclotron  irradia¬ 
tion  undergo  diffusion  with  reversible  trapping.  This  leads 
to  a  gradual  redistribution  of  the  deuterium  atoms  within 
the  metal  and  their  slow  evaporation  from  its  surface  -  an 
inverse  of  the  usual  adsorption  and  transport  mechanism. 

In  this  paper,  we  report  and  discuss  deuterium  and  oxygen 
atomic  distribution  profiles  followed  by  means  of  SIMS  and 
Auger  electron  spectroscopy. 

The  influence  of  hydrogen  on  fatigue  crack  growth  in 
structural  aluminum  alloys  is  discussed  widely  in  the  liter¬ 
ature^1-6^.  It  has  been  shown  that  humid  air  and  water 
vapor  increases  crack  propagation  rate  in  aluminum  alloys 

r 7  oi 

but  dry  hydrogen  has  no  effect  '  .  At  room  temperatures, 

diffusivity  of  hydrogen  in  metals  is  extremely  high  and  exceeds 

that  of  heavy  interstitials  like  oxygen  and  nitrogen  by  15  to 

[91 

20  orders  of  magnitude  .  Penetration  of  hydrogen  during 
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fatigue  crack  growth  could  be  deeper  than  the  predicted  depth 

for  simple  diffusion  using  the  sweep-in  concept .  The 

time  scale  of  hydrogen  presence  in  the  metal  during  fatigue 

is  certainly  much  longer  than  characteristic  time  for  hydrogen 

diffusion1  .  Therefore,  a  presence  of  trapping  sites  for 

hydrogen  must  be  expected  in  various  forms  such  as  Cottrell 

[9-141 

atmospheres,  interfaces,  vacancies,  voids,  etc.  .  It 

has  been  reported  that  hydrogen  implanted  into  the  single 
crystal  of  aluminum  diffuses  out  within  half  an  hour  , 
less  time  than  it  would  take  to  run  the  fatigue  experiment, 
thereby  confirming  the  need  for  trapping  sites  to  explain 
the  longer  presence  of  the  deuterium.  In  consequence,  a  con¬ 
centration  of  hydrogen  higher  than  the  solubility  level  in 
the  aluminum  alloy  associated  with  trapping  sites  could  be 
related  to  the  observed  increase  in  the  crack  propagation  rate 
during  the  fatigue  tests1  . 

Materials  and  Experimental  Approaches 

Specimens  of  2219-T87  and  7075-T6  aluminum  alloys  were 
prepared  in  four  different  ways.  The  preparation  prior  to 
deuterium  implantation  included  as  received  unrolled,  cold 
rolled  7-10%,  radiation  damaged  and  the  fracture  surfaces  of 
a  sample  fatigued  in  50%  RH  laboratory  air.  Surfaces  of  the 
as  received  unrolled,  cold  rolled  and  specimens  prepared  for 
the  radiation  damaged  received  light  polish  finish  using  a 
soft  cloth  polishing  wheel.  The  radiation  damaged  samples 
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were  subjected  to  the  neon  irradiation  with  a  beam  of  130  KeV 

15  2 

and  a  dose  8  x  10  ions/cm  at  a  projected  range  of  0.23  ym. 
Such  high  energy  ion  radiation  produces  point  defects  (inter¬ 
stitials,  vacancies)  and  vacancy  clusters'-'*'^.  Slightly 
higher  doses  of  impinging  neon  atoms  into  the  prepared  aluminum 
alloy  thin  disks  caused  overheating,  melting  and  destruction 
of  the  specimen.  The  specimens  were  cooled  to  room  temperature 
before  the  deuterium  implanting  process  to  be  described  below. 
Cold  rolling  produces  deformed  grains,  higher  dislocation 
density,  and  modified  precipitate  matrix  interfaces.  Thus, 
as  a  result  of  increased  dislocation  density  and  modified 
interfaces  a  cold  worked  material  has  a  higher  possibility  of 
trapping  hydrogen  than  an  unrolled  one.  The  fracture  surfaces 
from  the  fatigued  specimens  represent  a  graded  deformation 
process  with  the  largest  deformation  and  corresponding  trapping 
sites  expected  to  be  closest  to  the  surface. 

Due  to  the  high  background  level  of  hydrogen  in  the  AES/ 
SIMS  system  vacuum  chamber,  deuterium  as  an  isotope  was  used 

to  study  the  trapping  mechanism.  The  deuterium,  (D+)  implan- 

15  2 

tation  (beam  of  25  KeV,  and  a  dose  6  x  10  ions/cm  ,  at  a 
projected  range  of  0.24  ym)  was  done  in  an  AIM  210  Accelerator 
with  90°  magnetic  analyzer  and  energy  resolution  of  0.04%  FWHM. 
The  beam  was  perpendicular  to  the  implanted  surfaces.  The 
deuterium  implanted  specimens  were  transferred  into  the  AES- 
SIMS  vacuum  chamber.  Penetration  and  time  dependent  diffusion 
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profiles  of  implanted  ions  and  oxygen  were  obtained  using 
SIMS  and  Auger  Electron  Spectroscopy.  Due  to  the  experimen- 

_  -f- 

tally  observed  higher  yield  of  D  as  compared  to  D  ,  negative 
SIMS  was  utilized.  The  relative  yield  phenomenon  depends  on 
many  factors  such  as  ion  energy  as  was  observed  by  others  . 
In  addition,  platinum  masks  were  used  over  the  specimens  to 
minimize  the  contribution  to  the  SIMS  signal  from  non-uniformly 
sputtered  edges. 

Results  and  Discussion 

The  results  of  the  diffusion  of  the  ion  implanted 
deuterium  as  determined  by  inert  Ar  ion  sputtering  are  shown 
in  Figs.  1  and  2.  In  all  cases  the  concentration  after  24 
hours  is  greater  in  the  vicinity  of  the  surface  than  in  the 
interior.  However,  there  was  a  small  increase  in  con¬ 
centration  at  a  certain  depth  from  the  surface  for  neon 
predamaged  samples.  Assuming  an  average  inert  ion  sput¬ 
tering  rate  of  200A  per  min. ,  this  depth  closely  corresponds 
to  the  projected  range  of  both  implanted  ions.  This  increase 
in  concentration  of  deuterium  at  a  projected  range  is  due  to 
the  locally  induced  traps  formed  by  the  high  energy  neon  irra- 
dation.  It  is  highly  unlikely  that  near  surface  trapping 
sites  were  filled  starting  from  the  surface  during  the  implan¬ 
tation  since  no  such  attenuation  process  is  known  to  exist  for 
the  high  energy  incident  deuterium  ions.  The  binding  energies 
of  trapping  sites  with  deuterium  for  different  metals  are 
typically  five  orders  of  magnitude  smaller  than  the  kinetic 


'jfc  jj  iM 


-90- 


energy  of  impinging  ions ^ .  Moreover,  a  maximum  concen¬ 
tration  at  the  expected  penetration  depth  for  the  neon  predamaged 
specimens  is  well  defined  [Fig.  1,2].  Therefore,  the  initial 
implanted  distribution  is  that  expected  by  their  penetration 
to  the  projected  range  and  diffusion  of  deuterium  after 
implantation  must  be  responsible  for  the  ultimate  observed 
distribution.  The  relative  differences  in  the  final  concentra¬ 
tions  for  the  different  specimens  must  be  due  to  the  differences 
in  trapping  potential.  Some  quantitative  data  on  trapping 
energies  of  different  trapping  sites  is  available  in  litera- 
ture1  '  J .  Second  phase  particles,  solute  atoms,  low  angle 
grain  boundaries  and  microvoids  homogeneously  distributed  have 
been  considered  as  strong  traps  for  hydrogen  with  the  highest 
trapping  energies  on  the  order  of  -1  eV.  Microcracks,  large 
particles  and  high  angle  grain  boundaries  usually  heterogeneously 
distributed  are  considered  as  localized  traps  for  hydrogen 
having  an  increasing  damage  effect  on  the  material  and  trapping 
energies  lower  than  0.6  eV.  Dislocation  and  lattice  sites  are 
very  reversible  diffusion  sites  and  are  poor  traps  for  the 
hydrogen  with  a  trapping  energy  of  approximately  0.2  eV.  The 
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reported  trapping  energies  '  are  not  for  aluminum  alloys. 

In  the  case  of  aluminum  the  strongest  traps  should  be  voids  or 
microcracks  where  a  gas  can  form  which  would  not  tend  to  dis¬ 
sociate  on  aluminum  resulting  in  the  highest  trapping  energy . 
Fatigue  fractured  surfaces  are  expected  to  have  the  highest 
density  of  microvoids,  microcracks,  grain  boundaries  and 


dislocations  with  an  increasing  density  of  these  traps  towards 
the  surface.  The  accumulation  of  these  traps  at  the  surface 
is  expected  because  of  the  large  strain  field  formed  in  the 
vicinity  of  the  propagating  crack  due  to  the  crack  opening 
displacement  and  reverse  yielding.  In  case  of  cold  rolled 
specimens,  near  surface  accumulation  of  microcracks,  damaged 
boundaries,  precipitate  particles  and  high  dislocation  density 
are  known  to  be  sinks  for  hydrogen.  The  observed  lower  con¬ 
centration  of  deuterium  for  cold  rolled  specimens  in  comparison 
to  fatigue  fracture  is  not  surprising  due  to  the  presence  of 
lower  trapping  energy  sites.  The  as  received  unrolled  speci¬ 
mens  with  fewer  and  uniformly  distributed  trapping  sites  such 
as  second  phase  particles,  small  dislocation  density  and  very 
reversible  diffusion  lattice  sites  fall  in  the  third  category 
of  the  surfaces  capable  to  trap  implanted  deuterium.  The 
unrolled  specimens  and  the  neon  radiation  predamaged  specimens 
seem  to  behave  identically  in  terms  of  the  observed  near  sur¬ 
face  concentration  profile  for  deuterium.  The  light  mechanical 
polishing  of  the  unrolled  specimens  with  and  without  neon  pre¬ 
damage,  would  be  expected  to  introduce  limited  amounts  of  near 
surface  trapping  sites. 

To  explain  the  results  shown  in  Figs.  1  and  2,  sam¬ 
ples  can  be  separated  into  two  categories.  The  first  category 
includes  the  unrolled  specimens  with  and  without  neon  pre¬ 
damage.  In  this  case  the  source  of  the  observed  profile  is 


associated  with  the  oxide  surface  layer  acting  as  a  diffusion 
barrier.  The  second  case  is  for  the  fracture  surfaces  and  the 
cold  rolled  specimens.  In  these  cases  there  is  a  significant 
introduction  of  strong  traps  near  the  surface  but  extending 
interior  to  the  oxide  layer.  The  two  cases  are  discussed 
separately  below. 

In  the  first  approximation  we  can  consider  the  oxide 
layer  formed  at  the  surface  of  the  unrolled  and  neon  pre¬ 
damaged  specimens  as  a  barrier  for  the  outdiffusing  deuterium. 
The  initial  distribution  of  the  deuterium  just  after  the  implan¬ 
tation  can  be  represented  by: 


p (x,t  =  0)  =  lim 
t-*0 


/iDt 


exp 


(2l2iZ2oI 
'  4Dt 


-)  = 


<5 (x-xo) 


(1) 


The  increase  in  the  local  trapping  at  the  projected  range  in 
the  neon  predamaged  specimens  would  lower  the  total  flux  of 
the  diffusing  deuterium  towards  the  surface.  In  further  steps 
we  will  ignore  this  local  trapping  in  order  to  simplify 
general  form  of  the  solution. 

We  seek  solution  to  the  Fick's  equation 

2 

Mx,t)  =  -D  -2-4  (2) 


with  the  initial  condition  (1)  and  a  reflecting  barrier  at 
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3 P  (x,  t) 
Tx 


x=0 


0 


(3) 


It  is  customary  to  dispose  of  the  cumbersome  boundary  pro¬ 
blem  by  introducing  "mirror  image"  of  the  initial  distribu¬ 
tion.  This  yields  as  a  solution  the  sum  of  two  Gaussians: 


p(x,t)  = 


/rrDt 


(4) 


which  satisfies  boundary  condition  (3) .  In  the  region  signifi¬ 
cant  for  us  x>0,  p(x,t)  evolves  with  time  in  a  manner  shown 
in  Fig.  3.  The  solution  for  time  t4  is  what  we  observed  for 
near  surface  behavior  of  the  unrolled  and  neon  predamaged 
specimens  first  day  after  the  implantation.  The  second  day 
after  the  deuterium  implantation  no  traces  of  deuterium  were 
detected  in  unrolled  and  neon  predamaged  surfaces.  Therefore 
we  presume  that  near  surface  traps  in  these  specimens  must  be 
almost  uniformly  distributed,  reversible  and  very  weak  allowing 
fast  diffusion  of  implanted  deuterium  within  the  specimen. 
Furthermore,  the  thin  oxide  formed  at  these  surfaces  is  a 
strong  barrier  for  the  outdiffusing  deuterium. 

In  the  second  case  as  the  deuterium  diffuses  towards  the 
surface  in  the  cold  rolled  and  fatigued  specimens  it  gets  trapped 
in  the  high  binding  energy  traps  extending  their  presence  much 
deeper  than  just  the  oxide  layer.  Therefore  the  preceding 
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analysis  will  no  longer  predict  the  diffusion  profile  but 
a  more  complex  analysis  including  the  nonuniform  distribu¬ 
tion  of  traps  must  be  used. 

Traces  of  implanted  deuterium  were  observed  in  the  cold 
rolled  specimens  up  to  the  fifth  day  after  implantation  took 
place.  Afterwards,  up  to  the  14th  day,  only  fatigue  frac¬ 
tured  specimens  showed  remaining  concentrations  of  deuterium. 
Previously  discussed  strong  traps  present  in  the  surfaces  of 
these  two  specimens  are  responsible  for  such  a  long  listing 
presence  of  the  implanted  deuterium.  Deuterium  profiles 
obtained  with  SIMS  have  been  shown1  1  for  the  segment  of  spec¬ 
imen  fatigue  fractured  in  a  deuteriunated  water  vapor  environ¬ 
ment.  In  this  case,  dissociated  deuterium  entered  into  the 
bulk  of  the  material  through  the  free  metal  surface  exposed 
during  the  fatigue  process.  Deuterium  gets  trapped  in  a  num¬ 
ber  of  high  energy  trapping  sites  created  during  the  fatigue 
fracture  process.  The  immediate  formation  of  the  oxide  layer 
slows  down  the  back  transport  of  the  deuterium  to  the  surface 
where  evaporation  takes  place.  Deuterium  concentration  pro¬ 
files  obtained  from  these  specimens  and  fatigue  fractured 
implanted  specimens  showed  detailed  similarities.  Both  show 
a  high  concentration  of  deuterium  at  the  surfaces  and  a  signi¬ 
ficant  concentration  level  of  deuterium  extending  into  the 
material  much  deeper  than  the  observed  oxide  layers. 

This  observation  suggests  that  the  final  distribution  of 
deuterium  in  the  specimens  fatigue  fractured  in  the  D20 
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environment  and  ion  implanted  fatigue  fractured  surfaces 
strongly  depends  on  trapping  sites  formed  in  the  fatigue 
process. 

The  oxide-metal  interface  seems  to  play  another  impor¬ 
tant  role  for  deuterium  trapping  which  is  especially  pro¬ 
nounced  for  the  fatigue  fractured  specimens  of  the  7075-T6 
alloy.  Fig.  4  shows  the  correlation  between  oxygen  sputtering 
profile  obtained  using  AES  parallel  to  deuterium  profile 
utilized  from  SIMS.  At  the  oxide  metal  interface  a  signifi¬ 
cant  buildup  of  deuterium  ions  was  observed  suggesting  again 
an  oxide  to  be  a  barrier  and  a  oxide-metal  interface  serving 
as  a  strong  trap. 

Conclusions 

1.  Large  amounts  of  strong  deuterium  trapping  sites  close 

to  the  surface  are  created  during  the  deformation  associated 
with  fatigue  crack  propagation  process  in  aluminum  alloys. 

2.  Deuterium  profiles  of  D2O  vapor  fatigue  fractured  specimens 
were  closely  duplicated  by  trapping  inversion  experiment. 
Strong  dependence  of  deuterium  distribution  with  the  pres¬ 
ence  of  trapping  sites  created  during  the  fatigue:  fracture 
process  was  observed. 

3.  Surface  oxides  and  the  oxide-metal  interface  trap  deuterium 
slowing  down  the  transport  to  the  surface  where  evaporation 
takes  place. 
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4.  High  energy  neon  radiation  induced  damage  produces  weak 
trapping  sites  for  deuterium. 
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